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Aspects of the Turbulence Problem 
Survey Report 


By Hans W. LiEpMANN, Pasadena, California’) 


Second Part?) 


lV) GENERAL DISCUSSION 


1. Equations of Motion. Reynolds Equations 


We assume that the velocity components u,, the pressure #, the density @, 
the temperature T, etc. satisfy the general Navier-Stokes equations of motion, 
that is, we have 


Ogu; , OQUsUy, _ OTsK 
OF sna Ofae a 04,? (IV-1a) 
Oo F , 00S op 0 
Ot % on ee aaa [Tin Me + Gul » (IV-1b) 
00 , 90%, _ a 
a Ox, 0 (IV-1c) 
with 
. z OU, du; , 0 
Tk = —|p +4420) el bit Bae tae) (IV-14) 
oT 
Sr ee poe, TVv-1 
Ti OX ( e) 
JH=Z musth. (IV-1) 


A, wu are the two viscosity coefficients, « the heat conductivity, and h the 
enthalpy per unit mass of the fluid. 

A direct approach to the turbulence problem would consist in solving (LV-1) 
for a given set of boundary or initial values and to compute mean values over 
the ensemble of solutions. Even for the most restricted problem, turbulence 
of an incompressible fluid, this appears to be a hopeless undertaking due to 

‘the nonlinear terms in the equations. Thus, the standard procedure, following 
REYNOLDS’ classical studies, consists of averaging over the equations rather 
than over the solutions. In this way relations between averaged quantities are 
obtained. However, due to the nonlinear transport terms, the system of equa- 
tions obtainable in this way is undetermined since it is impossible to obtain 
an equation which involves only correlations of the same order, for example, 


double correlations. 


1) California Institute of Technology. 
2) First Part see ZAMP 3, 321 (1952). 
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Indeed, averaging (IV-la) for a field of flow consisting of a mean flow and 


turbulent fluctuations 
u,=U;,+u;,, %4=U; 


yields Reynolds equations which couple the mean flow U; to the quadratic 
mean values of the fluctuations. For example, for incompressible flow we have 
ete = (in — 0) (IV-2a) 


OX, Ox; 


00 U; 
ot 


a 


in which the ‘‘apparent’’ or Reynolds stresses now appear on the right hand 
side. Multiplying ([V-la) by w;, for example, and averaging we obtain an 
expression connecting the stress tensor «/ 1; with the expression containing 
the triple correlations w/w; u;, etc. The same procedure evidently applies to 
the general compressible case except that more terms and equations are 
involved. 

Hence, the method of averaging the equations can not lead to a final deter- 
mined solution since it evidently does not exhaust the content of the original 
equations. The averaged equations require additional relations which have to 
come from statistical or similarity considerations or from general physical 
considerations which allow a choice of the possible solutions to be made from 
the number of solutions compatible with the averaged equations. 


2. The Analogy with Kinetic Theory 


From Reynolds equations, for example, in the simple form (IV-2a), an 
immediate analogy between turbulence in a fluid and molecular motion can be 
drawn. The viscous stress tensor T;, expressed in terms of the molecular velo- 


cities c; is given by t,;, = —@ ¢; ¢, , thus (IV-2a) can be written 
00 U; do U; Op i) — esi 
Tres ~ OR, selon Ox, [0 6, C, + 0M; Uy]. (IV-2b) 


Ideas such as BoussINEsQ’s ‘exchange coefficient’, PRANDTL’s mixing length, 
etc. make use of the analogy between molecular and turbulent shearing stresses. 
It should be emphasized, however, that the computation of ¢;¢, from kinetic 
theory is reasonably simple only in the case of a perfect gas, that is, for weak 
molecular interaction. The interaction between turbulent elements, on the 
other hand, is not weak in all interesting cases. Indeed, if an analogy of turbu- 
lent mixing with molecular kinetics is drawn at all it has to be drawn with the 
liquid state rather than with the gaseous state. Furthermore, to consider the 
fluid incompressible implies that dissipated energy is small compared to the 
heat content of the fluid. That is, cc; <U,U, and the Mach number 
U?/c? ~ U*/c® is small compared to unity. For the turbulent fluid w/ u;, is not 
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necessarily small compared to U; U; and consequently the proper analogy 
refers to ‘‘compressible fluid flow” of this turbulent fluid. Indeed, the triple 
correlation terms of the turbulent fluctuations correspond to the heat flux 
tensor in the same way as the double correlation corresponds to the viscous 
stress tensor. For compressible flow the relation between heat flux and shear 
can not be neglected, in analogy to the relation between double and triple 
correlations in some problems in turbulence. 

One of the difficulties in dealing statistically with turbulence is evidently 
the fact that we have to deal with strongly interacting elements. This precludes, 
in most cases, the use of the general asymptotic laws of probability theory such 
as the central limit theorem which form the background to many of the results 
of statistical mechanics, for example, of perfect gases, black body radiation, 
etc. all of which are problems of independent or nearly independent systems. 

Furthermore, the fact that turbulence is a “‘secondary”’ structure super- 
imposed upon the molecular one means that we deal with dissipative systems 
throughout and not with a fluid in thermal equilibrium. The very simplest 
state conceivable is a stationary one in which energy is fed to the system at the 
same rate as it is removed from the system and transforr.ed into heat. 


3. Some Characteristics of the Equations 
Mathematical Models of Turbulence 


Turbulence is a phenomenon typical for large Reynolds numbers. It is thus 
important to ask what the general characteristics of the equation are for large 
Reynolds numbers. The typical phenomenon at large Reynolds numbers is the 
existence of boundary layers in a general sense, that is, of thin layers where 
most of the dissipation takes place. In this sense shock waves are included in 
the boundary layer phenomena. The reason for the existence of boundary layers 
is the fact that the stress terms in the equation of motion are the terms of 
highest order. These terms are multiplied by a factor inversely proportional to 
the Reynolds number. Physically this means the existence of regions of dimen- 
sion 6, say, within the fluid, in which the dissipation takes place and where 6 
tends toward zero as some power of », the kinematic viscosity, in such a fashion 
that the dissipated energy remains finite. This implies, for example, in the case 
of the shock wave 6 ~ », for the two-dimensional boundary layer on a solid 
wall 6~ y/2, in a plane jet 6 ~ y?/3, etc. (Similar considerations can be made 
for temperature layers.) The general tendency of the contribution of the quasi- 
linear terms in the equation of motion and the higher order stress term thus 
ts at high Reynolds numbers in the production of narrow 
zones of intense shear. The narrow shear zones caused by transversal waves 
of finite amplitude are vortex lines and vortex sheets. The random field of 
these transversal disturbances is called turbulence. In an incompressible fluid 


apparently consis 
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longitudinal waves do not exist. In a compressible fluid both types exist and if 
the Mach number is sufficiently high both the longitudinal and transversal 

phenomena, that is, sound and shock waves on one hand and turbulence on 

the other, are coupled, and energy is transferred between both fields. 

While shock waves exist in flow with one space coordinate, turbulence, at least 
in the case of incompressible flow and in the sense outlined above, is a strictly 
three-dimensional phenomenon. In two dimensions the balance between the 
regrouping or stretching of vortex filaments and dissipation is not possible. 
(Random fields made up of transversal shear waves are, of course, possible in 
_ two dimensions but they do not represent turbulence in the proper sense.) 

Expressed in terms of a Fourier representation of the velocity field the 
interplay of the nonlinear term with the stress term appears as a flux of energy 
from the lower frequencies to higher frequencies and the eventual dissipation 
at higher frequencies. In the course of time one effect of the nonlinear term 
consists in increasing the coefficients of the highest frequencies. The nonlinear 
terms thus contribute indirectly to the dissipation. 

BuRGERS [38]?) has extensively studied a mathematical model of turbulence. 
This model accounts for both the nonlinear term and the dissipation term. It 
is simplified mostly in the space dimensions since only one and two dimensions 
are considered. Hence, the model does not apply directly to turbulence. In 
the one-dimensional case it can be applied to a random arrangement of shock 
waves rather than to turbulence. It does, however, allow the study of the 
interplay between inertia and viscous terms. COLE [42], [43] and independently 
E. Hopr [49] have recently found a general solution of BURGER’s one-dimen- 
sional, nonlinear equation which facilitates an analytical discussion of the 
behavior of the solutions with respect to initial and boundary values. 

Recently an electrical model nearly corresponding to BURGER’s one-dimen- 
sional equation has been studied experimentally by BetTcHovy [37]. BETCHOV 
feeds noise with a white spectrum and Gaussian distribution into his circuit and 
observes the resulting output. The probability distribution of the response is 
skew and the joint probability of the output voltage and its derivative shows a 
striking resemblance to similar observations of turbulent fluctuations. 


4. Local Isotropy 


KoLMoGoROFF’s concept of local isotropy [51] can be considered as one 
of the most important general ideas in turbulence. 

The exchange of energy between eddies of various size or the flux of energy 
through a wave number space represents a cascade process. The energy is 
supplied essentially in the low wave number region and passed on and dissipated 
in the region of large wave numbers. The question arises: is it possible to make 


1) Numbers in brackets refer to the Bibliography, page 423. 
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an asymptotic statement concerning the form of the energy spectrum of the 

large wave numbers? This spectrum is due to a large number of random events, 
and we may hope therefore to be able to establish a general trend. The events 
are not independent and consequently the central limit theorem, or a similar 
statement intimately connected with independence, can not apply and it is to 
be expected that any asymptotic statement made will be much weaker than, 
for example, the central limit theorem. 

KOoLMoGorROFF introduces essentially three hypotheses; the first of which 
much broader than the other two. 

(1) The small scale motion is always isotropic. This is the concept of “local 
isotropy”. 

(2) The energy spectrum of the smail scale motion can not depend upon details 
of the flow; the energy spectrum E(k) or the space correlation function 
g(r) will depend upon the dissipated energy per unit volume and time « 
and the viscosity ». 

(3) If the Reynolds number of the problem is sufficiently large, the zone of 
dissipation and the zone of production of turbulent energy will be widely 
separated in the wave number space. Then there may exist a range of wave 
numbers which are in a state of local isotropic equilibrium but which are still 
not in the dissipation region. In this case E(k) will become independent of ». 
(2) and (3) imply a transfer mechanism which acts between neighboring 

wave number ranges. The contribution of a direct transfer of energy from 

one wave band into another far removed is assumed negligible. 

The simplest consequences of (1), (2) and (3) for the spectrum or correlation 
function are: due to (1) the pattern is isotropic, and one function suffices to 
describe correlation or spectrum; dimensional analysis can be applied to E(k). 
Characteristic length and velocity made up of v and ¢ are: 


P= yl g-Ut (IV-3a) 
c= pill gt (IV-3b) 


1, c can be interpreted as wave length and phase velocity of shear waves in a 
viscous fluid. (2) yields: 


Ekyecy et Clko e- 7"). (IV-4) 

For the subrange (3) » has to vanish from this expression and it follows that 
E(k) = const 7° R75! . (IV-5) 

This is the famous k~°? law, first given by KormocororF?) and independently 
1) Kormocororr, ONsAGeR, and von WeizsAcxer did not actually use the spectrum concept. 


They gave a result equivalent to ([V-6) for the correlation function, the result for the spectrum 
is actually due to HEISENBERG [83]. 
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discovered by ONSAGER [56] and von WEIZSACKER [66]). The corresponding 
double correlation function has the form 


g(r) = (1 — const e~ !8 7-2/8) w2, (IV-6) 


Clearly, (2) implies that all quantities associated with eddies within the local 
isotropic zone are expressible in terms of (IV-3) and that they are similar. 

Kotmocororr’s results have an asymptotic nature. Hence for a comparison 
of, say, (IV-5) or (IV-6) with experiments the range of Reynolds numbers in 
which the asymptotic laws can be applied has to be established. Only after this 
has been done can the adequacy of the general ideals underlying the hypothesis 
be checked with experimental results. The status of an experimental verifi- 
cation appears to be as follows: 

(i) Local Isotropy: The existence of local isotropy in shear flow appears 
rather well established. The most convincing results are the measurements of 
CorrsIN [46] and of LAUFER [52] in a jet and channel respectively. CORRSIN 
and LAvuFER show that the contribution of small scale eddies to the apparent 
shear vanishes faster than their energy. Extensive results concerning the dissi- 
pation, skewness, etc. are due to TOWNSEND [61], [62]. 

(11) Similarity in the Viscous Region: Results of STEWART and TOWNSEND 
[99] in isotropic turbulence show that the high frequency components of the 
spectrum are similar in the sense of (IV-4). 

(iii) Nonviscous Subrange: Measurements of the spectrum or correlation 
function in isotropic turbulence behind grids have been carried out to establish 
whether or not a subrange of the form (IV-5) or (I[V-6) exists within the range 
of wind tunnel experiments. At the Reynolds numbers obtained so far in wind 
tunnels the results show that no extensive subrange of this form exists [61], 
[91], [99]. This has been most clearly demonstrated by STEWART and Towns- 
END, who estimate that the Reynolds number based on the grid mesh required 
to obtain a nonviscous subrange should be of the order of three million or 
roughly ten times larger than hitherto investigated. 

Hence, wind tunnel measurements have not yet lead to a decision concern- 
ing KoLMoGoroFr’s nonviscous subrange and the form of the corresponding 
spectrum and correlation function. 

Recent investigations by MacCreapy [54] of the spectrum of atmospheric 
turbulence near the ground give support to the k~ °/* law. MACCREADY’s measure- 
ments indicate that local similarity exists for wave lengths up to 160 cm. 
Reasonable agreement with the k~°/* law or the corresponding 7~2/? law was 
found for even larger eddy dimensions up to an order of a hundred meters. 
While the results of MACCREADY are very suggestive, the accuracy of measure- 
ment for atmospheric turbulence is naturally less than for wind tunnel investi- 
gations and hence MacCreEapy’s result can not yet be taken as entirely con- 
clusive evidence. 
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Some astrophysical results concerning the existence of the k~*/? law have 
been given by vON WEIZSACKER [67] and von Horner [48], but these results 
are rather inconclusive due to the difficulty of proper interpretation of the 
observations. 

The concept of local isotropy by its very nature does not contribute 
directly to an understanding of the transport phenomena. The Kolmogoroff 
region contains little energy but, at high Reynolds numbers the whole dissipation 
_ should occur within the region and hence the dissipation terms are identical with 
the simple expressions for isotropic turbulence. However, there are physical 
situations in which one is mainly interested in the small scale motion even if it 
does contain comparatively little of the total turbulent energy. Typical examples 
are the gust problem for the motion of an airplane in atmospheric turbulence, 
scattering of waves due to turbulence, etc. For problems of this type where 
the general turbulent motion is either unknown or not even interesting, a general 
asymptotic statement is of paramount importance. 

BATCHELOR’S work on axisymmetrical and general nonisotropic but homo- 
geneous turbulence should be mentioned here [36]. BATCHELOR shows that the 
pressure-velocity correlation terms in the equation transfer energy at one 
wave number from one velocity component to another and thus contribute 
a trend toward isotropy. In isotropic turbulence these terms vanish. 


5. Turbulent Flows with Secondary Structure 


In recent years the importance of the existence of a secondary, large scale 
structure in turbulent shear flow has become apparent. CorRsIN [44] and 
TowNnsEND [60] found that the flow near the outer edge of a jet or wake is only 
intermittently turbulent. TowNsSEND has shown the importance of this motion, 
which is of a scale comparable to the width of the wake, upon the momentum 
and heat transfer in the wake. Recent experiments by TowNsENP [65] and 
KLEBANOFF and Dieu [50] demonstrated the same intermittency near the 
edge of a boundary layer. Photographs of the wake of high-speed projectiles 
show that the turbulent wake has a large scale structure superimposed upon the 
fine scale turbulence quite similar to the wake structure at low speeds. 

Intermittency and thus the existence of elements of a very large scale seem 
to be typical for turbulent flows with free boundaries. 

While the large eddies ordinarily found in intermittently turbulent flow 
appear statistically distributed, there do exist a number of cases in which a 
regular or nearly regular motion of large scale superimposed upon turbulent 
flow has been observed. Par [57] and Mac Puatt [55] found that the three- 
dimensional vortices which initiate the laminar instability in the flow between 
rotating cylinders (TAYLOR [100]) persist even if the flow has become fully 
turbulent. In recent measurements on a vortex street behind a cylinder RosHKo 
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[58] found a similar result: above a cylinder Reynolds number of about 150 
the wake flow is essentially turbulent with superimposed, nearly equally spaced 
vortices. It is very likely that similar patterns exist in other situations, for 
example, in zones of large thermal instability. 

The importance of these large eddies lies in their bulk transfer of turbulent 
fluid and their consequent influence on transfer properties and on the energy 
balance. Due to the size of these elements, which is comparable to the physical 
dimension of the flow, for example, the width of the wake, it is not possible to 
account for their behavior in terms of local quantities. For the case of the wake 
TOWNSEND [63] has discussed the possibility of a quasi-equilibrium taking into 
account this secondary structure. While there may be doubt about the details 
of this secondary structure there is no doubt that the large scale motion can 
not be ignored for a great many problems including problems in sound pro- 
duction from jets, combustion, etc. 

If we consider turbulence as a sort of secondary structure superimposed 
upon the flow of an ordinary viscous fluid we can look at the large scale motion 
again as a superstructure on the turbulent motion. The interesting and signi- 
ficant feature of this hierarchy is the fact that the length scales are strictly 
separated. This feature is probably intimately connected with the form of the 
equation of motion and again related to the boundary layer phenomenon. 


V. ISOTROPIC TURBULENGE 


In dealing with turbulence in nonuniform mean flow we are faced with 
a typical transport problem. The intensity of the turbulence varies in space 
and possibly in time and the gradient of the mean velocity introduces a preferred 
direction into the problem. We have to deal with a nonhomogeneous and non- 
isotropic state of turbulence. The most fruitful results in turbulence resulted 
from G. I. TAyLor’s introduction of the concept of isotropic and homogeneous 
fields of turbulence. Consider, for example, a box filled with a fluid. At the 
time ¢= 0 the fluid is stirred up in a suitable way, say, by dropping a grid © 
through the box. The turbulence created in this fashion can be assumed to 
become isotropic and homogeneous in space, but it decays, so that the field is 
not stationary. It is possible, of course, to conceive a setup such that energy 
is continuously fed into the system by the stirring motion and removed again 
in the form of heat created by viscous dissipation. In considering this latter 
problem both the similarity and the essential difference between a box “filled 
with isotropic, homogeneous turbulence” and a gas or radiation in thermal 
equilibrium become evident. The similarity consists in the lack of over-all 
transport of momentum, heat, etc.; only the interactions between the elements 
making up the fluid are important. The difference is that in the turbulent 
system there is an energy flux which must be maintained by an external agent, 
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that is, the system is dissipative. To be sure, we may well question the possi- 
bility of such a setup. For example, the initial motion set up by a stirring 
mechanism can not be isotropic and homogeneous. Hence, we assume a trend 
toward isotropy rapid enough that between the production of turbulence and 
the dissipation into heat there exists a sufficiently long time to set up an 
isotropic homogeneous state. The reasoning here is evidently based upon ideas 
similar to those used in the concept of local isotropy. It may be expected that 
an isotropic state is reached more rapidly the more uniform the stirring 
mechanism. 

Granted that a state of isotropic and homogeneous turbulence is possible, 
we can set up the kinematics of the field and the equations interrelating the 
mean values. The general scheme to introduce correlation tensors and to set 
up the equations from there is due to voN KARMAN and Howartu [85]. In the 
following we will restrict the discussion to the incompressible case. 


1. Correlation Functions 


The field of isotropic, homogeneous turbulence is described by the two-point 
correlation tensors. Take two points P,(R) and P,(R-+ 7) in the field and form 
the tensors 


Pressure Correlation 
n(r) = 2(R) (R+7) =P’ (V-1a) 
Velocity, Double Correlation 


Rjj(7) = ui(R) u(R +7) = 05 05 (V-1b) 
Velocity, Triple Correlation 
Tegel?) = a(R) uj(R) (R$ 7) = me 15 tH « (V-10) 
For isotropic and homogeneous fields 
m=n(r), Rin= Rr) 17x + Rol) Ore» | (v-2)3) 
Tyin = Tilt) 12761 + Tar) 12 85% + Ta) 15 Ocn + Ta”) 7x 5:3, ete. | 
For an incompressible fluid the continuity equation gives 
Ous 
Oa » 
ei oR OT, 
p(R) u,(R+7) =0; Da 0; ov, Bee, Os (V-3) 


1) f(r), h(r), etc. is sometimes written for f(r, #), h(r, t), etc. Similarily t’, h’, etc. denote deriva- 


tives with respect to 7. 
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Hence R,,;, T;;, can each be expressed by a single scalar function of 7. Following 
the Kérman-Howarth notation we use the following functions: 


Ruf, 0, O}= f(r); R,{ 0,7, O}= g(r), ae 
Tur{7,0,0}=kr); Thf 0, 7, O}=hA(r); Tro{7, 0, 0}= g(r) . 
From continuity 
g=ft$f; k=2h; g=—(h+ 5h). (V-5) 
It is also useful to introduce the trace of R;, 
Ri = Ra = F2¢e (V-6) 
and 
Tina TO: TW=5- (V-7) 


Furthermore, the products of the derivatives of the velocity components can 
be expressed as derivatives of the correlation tensors. This is shown in exactly 
the same way as in the case of one-dimensional stochastic processes (II-11) 


= (Sa 8 du;\ (0? Ri; ) 
OX) OS en a ~~ es O%m /r=0 


' (V-8) 
= #"(0) [2 d1m 815 — (8:2 94m + dim 950] 
which leads to the expression for the dissipation function ®: 
£ OMS Nn " 
Na bap soe —15 wf" (0) =—5 yu R"(0). (V-9) 
The mean square vorticity Q? is also easily found from (V-8) 
= -15 nj") = 2 (V-10) 


- 
2. Spectrum Functions 


Using the same formalism as in Section II we can introduce the Fourier 
transforms of the correlations functions, the spectra. This was first done by 
TAYLOR [103] for the so-called ‘one-dimensional spectrum”, that is, the 
transform of /(7). HEISENBERG [83] introduced the three-dimensional spectrum 
and KAMPE DE FERIET [15] and BATCHELOR [36] made use of a general spec- 
trum tensor ®;; obtained by a Fourier transformation of R; ;. Thus 


@,,(k) = 


(V-11) 
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where k denotes the wave number vector. For isotropic turbulence, @®,; must 
be of the form 


D; j(h) = D,(R) he hy + Dal) 5. (V-12) 


Continuity for an incompressible fluid means the existence of transversal 
waves only, hence it follows that 


@;;k; = (Ok? + ®,) k; = 0 
and thus 


©,; = ©, k? |=ir Sais (V-13) 


from which HEISENBERG’S spectrum E(k) is obtained by 


E(k) =22k? ©,,= —82h'@,. (V-14) 
Thus 
= E(k hy ky 
@5(h) = he [8,5 — HR). (V-15) 


F(k,), TAYLOR’s one-dimensional spectrum, results from this by noting from 
(V-14), (V-11) and (V-7) that 


= Decal eG ry) krsinkr dr (V-16) 


thus 
E(k) => ~ [3 F"(hy) — hy F' (Ay) |p, 2 - (WeN7} 


The relation between TAYLOR’s and HEISENBERG’S form of the spectrum is 
quite analogous to the relation between, say, the Maxwell distribution function 
for one component of the velocity g(u,) and the distribution function for the 
absolute velocity p(c). The latter sorts the molecules with respect to their 
kinetic energies while the former sorts only according to the contribution of 
one component to the kinetic energy. Thus (0) + 0 while vo) = 0 and 
similarly F(o) + 0 while E(o) = 0. y(c) and E(k) have more direct physical 
significance, however y(u,) and F(k,) are more easily accessible to measurement. 


3. The Equations of Motion 


The Karmdan-Howarth equation establishes the relation between R;,; and 
T,;, and their derivatives. The equation can be written in the form 


of 2 Oh v4 2» 0 ( 4 -) (V-18 
Sieiie Fier roy" or}? : ) 
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where /(7, 4) and h(r, ¢) are the correlation functions defined in (V-5). Another 
useful form is 
OR S42 feat bare ase a aa 
oe ee T) = OF (r? R), ( ) 
where R(r, ¢) and T(r, t) have been previously defined. Multiplying by kr sink r 
and integrating over all ry we have, with (V-16) for the function E(A, ?) 


nef Ory ey ara 2 BR. (V-20) 


4. Invariants 


From (V-18) and (V-19) follows an invariant given first by LoITsIANSKII 
[93]. Namely, multiplying the equations by 74 and 7? respectively and inte- 
grating over all space yields 


if r* f(r, t) dr=0, (V-21) 
ss " r® R(r,t) dr =0 (V-22) 


provided that /, 7, k, T vanish sufficiently fast for large values of r. Further- 
more, since R = (r3 7)’/r3 


ter R(r, t) dr=0. 
The integral in (V-21) 


/ f(r) rAdr=A (V-23) 


is called LoITsIANsku’s invariant. The existence of such an invariant is a simple 
consequence of the equation of motion. The continuity, momentum and energy. 
equations in fluid dynamics each have the form of a continuity equation 
relating the time derivative of a function to a divergence term. The divergence 
terms can be written in the form 


ay; (OW) mn] = Fe O77, =F @ rs), ete. 


and it becomes evident that for homogeneous fields invariants will exist, For 


the velocity field this is Lorrstansku1’s invariant, for the temperature field in 
isotropic turbulence CoRRSIN [79] has given the appropriate invariant 


ae if POX dr= (V-24) 
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where # is the temperature fluctuation. For a compressible fluid a similar result 
holds for the density fluctuation s: 


apf 7259 dr =0 (V-25) 


as given by CHANDRASEKHAR [77]. 

The consequence of the existence of the invariant A on the spectrum F(R) 
follows directly from (V-16) and (V-20). A, being a fourth moment of f, is 
essentially determined by the largest elements of turbulence and hence it must 
determine E(k) in the neighborhood of k = 0. Developing (V-16) and (V-20) 


into powers of & yields 
BIS Al 


since 
1 , 
Ri) = 4s 8A. 
Similarily, remembering that T(r) = —(r4 K)'/r? (V-20) yields 
dE ke : 
ces / rT (r,t) dy -2vk®E (for small &) 


consistent with the statement 0A/0t = 0. 

Consequently, the spectrum function F is proportional to k* near the origin 
and the coefficient of k4 is time independent. This result was found by LIn [92] 
and BATCHELOR [36]'). 

We again draw an analogy to the Maxwell distribution of velocities or the 
Plank distribution of radiation. In both cases we can also obviously express the 
distribution functions in such a form that their behavior near the origin 
remains invariant as the gas or the cavity cools down with time. However, 
if this is done the whole distribution remains similar, which is typical for a 
linear system and quite different from the turbulent case in general. 


5. The General Problems in Isotropic Turbulence 


The discussion of invariants and of the behavior of E(k) for small k, the k* 
law, essentially exhausts the results which can be obtained from the Kaérman- 
Howarth and related equations without the introduction of additional prin- 
ciples or assumptions. The general problem posed to a theory of isotropic 
turbulence is the determination of the complete functions E(k, ¢), f(7, 4) or 
D, i(k, /), etc. It is immediately questionable whether such a broad problem 
makes sense. The necessity of a stirring mechanism has to be considered some- 
how and it is by no means evident that turbulence created in some manner, 


1) BarcueEor shows that the k* law requires homogeneity only. The coefficient of k4 is in this 
case still invariant in time but is a tensor term and not a simple scalar like in the isotropic case. 
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for example, in a box, will tend toward a unique function E(k, ¢). Indeed, the 
existence of the Loitsianskii invariant does establish a connection with the 
stirring mechanism for all times. For example, the anisotropy introduced by 
the stirring mechanism can be recovered in the latest stages of decay, as shown 
by BATCHELOR und Stewart [70]. Hence, while it may be futile to ask for the 
complete spectrum or correlation function without a complete description of 
the stirring mechanism, it makes sense to ask for any possible general form the 
spectrum may have in wave lengths far removed from the ones excited by the 
stirring mechanism and for times long with respect to a relaxation time. It may 
be sufficient in these cases to describe the stirring mechanism by a character- 
istic length only, for example, the mesh size in the case of a grid. 
The approaches taken to a solution of the problem can be classified into 
groups. 
(a) Similarity considerations and dimensional analysis. 
(b) Attempts to express the triple correlation or the corresponding spectral 
transfer integrals in terms of the double correlation or the spectrum E(k). 
(c) Discussion of limiting cases. 


6. The Final Stage of Decay 


' The difficulty in discussing the Karman-Howarth or related equations 
arises naturally from the transport term. One may thus first consider a case 
where this term is negligible, as has already been done by voN KARMAN and 
HowartH [85] and MILLIONSHCHIKOV [94] and more completely by BATCHELOR 
and TowNSsEND [72]. If the turbulence has decayed far enough for the transport 
term to be negligible (V-18) or (V-20) may be integrated. The appropriate 
solutions are 


E(k, t) ==.constk4e-*”**, (V-27) 
f(r, t) = (82 vy 98 [Hr — §, ty pe 8"! 54 ds (V-28) 
and thus 
Vee LA #)-5/2 
u Win (wt) 77". (ast->oo) (V-29) 


(V-27) and (V-29) are in good agreement with experiments [72], [91]. 


7. Self-Preservation 


In the Kolmogoroff region the spectrum and correlation functions depend 
upon time during decay only indirectly through the dissipation e. For the larger 
part of the energy spectrum that contains much of the turbulent energy, time 
has to enter explicitly. Forms of the spectrum or correlation function which 
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depend upon time but preserve their shape are called self-preserving. Von 
KARMAN and LIn [87] and Stewart and TowNsEND [99] have recently discussed 
problems of self-preservation in detail. HEISENBERG [82] used the concept of a 
quasi-equilibrium using his form of the exchange term (see below). Von KARMAN 
and Lin put E(k, t) = C k* ©(k/k,) where k, is a function of t. This expression 
leads to a decay law of the form 


U2 my £1017 | (V-30) 


HEISENBERG puts E(k, t) = E)(t) E(k/k,) consistent with an initial decay law 
of the form 


u* we tI, (V-31) 


While von KArMAN and Lin’s expression for E satisfies the proper condition 
E ~k* for k +0, HEISENBERG’S expression yields E ~k for small k. The 
t-197 law proposed also by KotmocororF [88] and FRENKIEL [80] has not 
much experimental support while the ¢-1 law first discussed by BATCHELOR 
and TOWNSEND appears to fit the experimental results in the initial stage of 
decay remarkably well. 

For a detailed discussion reference is made to the papers by voN KARMAN 
and Lin, and by STEWARD and TOWNSEND. 


8. Form of the Exchange Term 


The spectral equation (V-20) integrated from 0 to # reads 


k 


k k oo 
0 Bae: i Oye sink 7 me [Re 
ef Bb) dk= =| (rey DEY dr dk ag R°E(k,t)dk. — (V-32) 
0 Oyo) 


The time rate of change of the energy contained in the range of wave numbers 
up to k is due to exchange with elements outside the range, as expressed in the 
second term, and to viscous dissipation. To make the equation determinate the 
exchange term has to be expressed in terms of E(k, ¢) and k. None of the forms 
of the exchange terms proposed so far have been completely successful. HEISEN- 
BERG’Ss form, or a more general similar expression due to voN KARMAN, seems 
to have been the most succesful so far. HEISENBERG interprets nonviscous loss 
of energy from the range of waves numbers 0 to # as due to a turbulent viscosity 
in the sense of BOUSSINESQ’S exchange coefficient. Using dimensional reasoning 
he then puts 

2 e sink y ff oe NArz i A V-33) 
of fiery pane edhe sh) (<x) ak | k* E dk. (V-33) 


0 0 0 
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This form of the exchange term yields the k~°/* law in the nonviscous subrange 
of Ko_mocororr, as any dimensionally correct forms do, but leads to a law 
of the form k-7 in the viscous subrange of KormocororF. This latter result is 
at variance with experiments ([106] and [91]) and also a priori unlikely since 
any power law E(k) ~ k-” for large k leads to infinite moments of E of order 
n — 1. This in turn means that the velocity components have no mean deriva- 
tive of order (n — 1)/2, which is not likely. PRouDMAN has recently discussed 
the consequences of HEISENBERG’S assumptions [96]. 

The proposed forms of OBouKHOV [95] and KovaszNay [89] are less suc- 
cessful than HEISENBERG’S. 


9. Temperature Fluctuations in Isotropic Turbulence 


CorRSIN [79] has studied the case of temperature fluctuations in a field of 
isotropic, homogeneous turbulence of an incompressible fluid. In this case he 
obtains an equation similar to the Karman-Howarth equation obtained for the 
temperature correlation function 


a 


O(R) O(R +7) = 6(r) (V-34) 


and the temperature velocity correlations 


oe 


O(R) O(R + 7) wi(R + 7) = O,(r) . (V-35) 


Similarily equations for the spectrum can be written. The difference due to the 
scalar character of # as compared to the vector character of «; becomes apparent 
in the form of the spectrum for small k. The temperature spectrum begins with 
k* corresponding to the invariance of 


| ‘OW rdr (V-36) 


during decay. The discussion of the final state of decay, the subranges, etc. is 
similar to that for the velocity field, as shown in CorRsIN’s paper. 


10. Pressure Fluctuations in Isotropic Turbulence 


The pressure of an incompressible fluid / satisfies the equation 


° 
Reh 
Ox; Ox; 


Ill 


pop i. (V-37) 
(V-37) is a very simple form of the general type of (III-6). Thus we immediately 
obtain from (III-11) the relation between the spectra of the pressure and the 
second velocity derivatives. Denote by J7(k) the pressure spectrum and by F(k) 
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the spectrum of f in (V-37). Then 


il 


IT(t) = = F(R) (V-38) 


| since ["(k) = k~?. (V-38) has been given by BATCHELOR [75]. BATCHELOR intro- 
duces the spectrum tensor ®; ;; ,,() corresponding to the quadruple correlation?) 


Pisim(?) = u4(R) u;(R) (RK +7) thy (R + *) — u,(R) uj(R) + u,(R+ 1) bm(R + 7), 
(V-38) then becomes 


> 


RA TT(R) = be hy by hin Desrm(®) (V-39) 

and 
p? =42 02 ih R2 IT (k) dk, (V-40a) 
(grad p)? = 4.2 0? | HTe) dk. (V-40b) 


However, to evaluate p; 51m or ®;; 7 it is necessary to introduce further assump- 
tions. BATCHELOR and also earlier HEISENBERG [83] assume the relation between 
the quadruple and double correlations to be the same as the one resulting from 
a normal joint probability of u(R) and u(R + 7). With this assumption JZ(k) 
can be related to the energy spectrum function E(k) of isotropic turbulence and 
approximately evaluated. 

To relate the statistical properties of the substantial acceleration in hydro- 
dynamics to the derivatives at a fixed position we must know the properties 
of the pressure gradient. A study of the pressure field is thus very important 
for the relation between the Eulerian and Lagrangian correlations and thus 
for the diffusion problem in turbulence. Besides BATCHELOR’s paper [75] 
reference is made to recent work of CoLis [78] and UBERoI and CorrsIN [107]. 
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Zusammenfassung 


Die Arbeit gibt einen Uberblick iiber einige Probleme der neueren Turbulenz- 
forschung. Nach einer allgemeinen Einleitung werden zuerst elementare, fiir die 
Turbulenztheorie wichtige Resultate der Theorie homogener, stochastischer Vor- 
gange diskutiert. Im dritten Teil werden sodann Beispiele linearer Systeme be- 
handelt, bei denen die dusseren Krafte von turbulenten Schwankungen herriihren. 
Im vierten Teil werden turbulente Transportphanomene gestreift, und im fiinften 
Teil wird ein kurzer Uberblick iiber Resultate und Fragen der Theorie der iso- 
tropen Turbulenz gegeben., 
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Zwei Integralgleichungen erster Art, die sich mit Hilfe 
Mathieuscher Funktionen ldsen lassen’) 


Von JOHANNES DORR, Darmstadt?) 


1. Aufgabenstellung 


Wir betrachten die folgenden beiden Integralgleichungen erster Art: 


g(x) = | H() Zolb |v — 9) dy (1 
und yp 
uaa helt iba iy 
w(x) = f p(y) Za Z, (ble — y|) dy 2) 
mit : 
Z,(et) = A Inlet) + BN,(a) 3) 


Darin bedeuten die J, bzw. N,, die Besselschen bzw. Neumannschen Funktionen 
n-ter Ordnung, und A und B sind beliebige komplexe Konstanten. 

Die erste Integralgleichung tritt auf bei der Untersuchung von Beugungs- 
erscheinungen von Schall- oder Lichtwellen am geraden Spalt, und zwar muss 
dann Z,= H2 gesetzt werden, wenn man die Zeitabhangigkeit durch 
exp (+7 #) ausdriickt. Die zweite Integralgleichung ist mit 7, = Hj? der we- 
sentliche Bestandteil der Integralgleichung fiir den Druck # und den Abwind w 
des harmonisch schwingenden diinnen Profils im kompressiblen Medium. Wird 
die Zeitabhingigkeit durch exp (—? @ /) ausgedrtickt, so sind die Hankel- 
Funktionen zweiter Art durch die erster Art zu ersetzen. 

Die beiden Integralgleichungen (1) und (2) sind eng miteinander verwandt. 
Die zweite geht aus der ersten durch Differentiation nach x hervor. Die Losung 
der zweiten Integralgleichung lasst sich deshalb relativ einfach aus der Losung 
der ersten herleiten. Es wird daher vorerst nur die erste Integralgleichung 


untersucht. 


1) Ausfiihrliche Wiedergabe eines auf der GaMM.-Tagung (Freiburg i. Br. vom 29, bis 31. Marz 


1951) gehaltenen Vortrages. 
2) Institut fiir praktische Mathematik der Technischen Hochschule Darmstadt, 
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2. Die Eigenfunktionen der Integralgleichung (1) 


Mit der Substitution 
x=-—cosB, y= —cosa (4a) 


nimmt (1) die folgende Form an: 


G() = | He) sina Z, (k |cosa — cos f|) da . (4b) 


Setzt man noch 
H (a) sina = F(a) , 
so wird daraus 


> 


G(B) = | Fe Z (k |cosa — cos B|) da . (4c) 


Damit hat unsere Integralgleichung eine Form erhalten, die den Anschluss an 
bekannte Ergebnisse erlaubt. WHITTAKER?) hat die Giiltigkeit der Beziehungen 


+7 


cat aVadae fe ceq(ax) Jg[k (cosa — cosf)] dx (5) 


bewiesen, wobei die ce,(f) die geraden Mathieuschen Funktionen sind, die der 
Differentialgleichung 


G"(B) + (q — k® cos® B) G(B) = 0 (6) 


gentigen. Man gewinnt die Beziehungen (5), wenn man von den ebenfalls von 
WHITTAKER gefundenen Beziehungen 

Ce, (B) = py / ce, (x) exp (t k cosa cosB) dx (7) 
ausgeht, fiir ce,(«) wieder die Integraldarstellung der Form (7) einsetzt und 
dann die Reihenfolge der Integrationen vertauscht. 


Es wird nun gezeigt werden, dass nicht nur die Relationen (5) , sondern 
auch noch die Relationen 


ce,(B) = a2 if cen(x) Ny (& |cosa — cos B|) dx (8) 


fiir alle 7 = 0, 1, 2, 3, ... giiltig sind. Damit ist zugleich gezeigt, dass auch die 
allgemeinere Integralgleichung (4c) die geraden Mathieuschen Funktionen als 
Eigenfunktionen besitzt. 


1) E. T, WuitTaKer und G.N. Watson, A Course of Modern A nalysis, 4. Aufl. (University Press, 
Cambridge 1927), S. 426. 
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Der Versuch, die Beziehungen (8) ebenso wie die Beziehungen (5) aus 
andern bereits bekannten Integralgleichungen der Mathieuschen Funktionen 
herzuleiten, gelang nicht. Wir nehmen deshalb funktionentheoretische Metho- 
_ den zu Hilfe und betrachten dazu den Ausdruck 


vB) = | 9) Kl. B) ax (9a) 
mit a“ 
K(a, B) = NM (& | cosa — cos§)) - (9b) 


Wir wollen voraussetzen, dass sich y(«) als Funktion der komplexen Variablen 

a auffassen lasse und in dem Streifen —z < Re [a] < +a regular analytisch 
sei. Der Kern K(«, 8) ist aber nicht ohne weiteres in diesem Streifen der 
a-Ebene als regular analytische Funktion deutbar. Wir betrachten deshalb die 
Funktion 


Ka, B) = Ny [k (cosa — cos£)| | (9c) 
und fassen hier « als komplexe Variable auf. In dem Streifen 
—z SRelal Sin (10) 


hat S(a, 8) fiir reelle B-Werte (~ x <8 S+2) auf der reellen Achse der 
a-Ebene zwei logarithmische Verzweigungspunkte an den Stellen «, = + f. 
Wir machen S(a, 6) in dem Streifen (10) der «-Ebene eindeutig, indem wir 

von dem Verzweigungspunkt «,= + |f| einen Schnitt entlang der reellen 
Achse nach rechts und von dem Verzweigungspunkt « = — || einen Schnitt 
entlang der reellen Achse nach links fiihren. Dann kann man auf Grund der 
Umlaufrelation fiir den Logarithmus an Stelle von (9a) auch schreiben: 


5 [vo 36 Bias | 0 R(x, B) de (11) 


eh 


Dabei bedeutet das Zeichen ~ bzw. v, dass der Integrationsweg oberhalb 
bzw. unterhalb der Verzweigungsschnitte von — a nach + a zu fiihren ist. 
Wir fordern vorerst 
Sarah ee (12) 


Dann ist (a, B) tiberall auf beiden Integrationswegen regular analytisch und 
gentigt, wie man leicht priifen kann, folgender partieller Lessee WA At 


oe ta #2 cost) R= 2S + (q— BY cosa) B- (13) 
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Darin bedeutet q¢ einen vorerst beliebigen Parameter. Da unter der Bedingung 
(12) die Integranden von (11) auf beiden Integrationswegen regular analytisch 
sind, darf man, wenn man den Ausdruck (11) nach f differenziert, die Reihen- 
folge von Integration und Differentiation vertauschen. Auf Grund von (13) 
findet man deshalb aus (11) nach zweimaliger Anwendung der partiellen 
Integration: 


yp" (B) + (¢ — k? cos? B) p(B) 
= 5 | [o"@) + @— cost a) o(a)] Ka, B) de 


+> | [p"(a«) + (¢ — k® cos? x) p(a)] RK(a, B) da 


a= +H 


+ [p(@) so K(@, B) — 9'@) K(a, B)| 


Die beiden Integrale lassen sich zusammenfassen. Ausserdem ist 


4 Kl, p)=0 fir a=4+x, K(—2,f)=K(+a). 


Cm —I7 


Damit erhalt man schliesslich 


yp" (B) + (q — k? cos? B) p(B) 


+7 


= | [e"(2) + q— #8 cos* x) y(a)] K (a, p) da (15) 


+ K(x, B) [p'(— 2) — 9'(+a)]. 


Aus (9) folgt, da der Kern K(x, 8) eine gerade Funktion in « ist, dass fiir eine 
ungerade Funktion g(«) die zugeordnete Funktion y(f) identisch Null ist. Des- 
halb wollen wir fiir g(«) nur gerade Funktionen zulassen. Es sei also: 


P(x) = y(—a) . (16a) 
Wir fordern nun ausserdem 
y(-2) =9'(+2), (16b) 
p" (%) + (¢ — k? cos? «) p(x) = 0. (16c) 


Durch diese drei Forderungen werden die geraden Mathieuschen Funktionen 
erster Art ce,(«) definiert, deren Existenz aber an die Bedingung gekniipft 
ist, dass bei vorgegebenem k-Wert der Parameter q spezielle Werte 


7 = Ink) (16d) 
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annimmt. Wenn alle vier Forderungen (16) erfiillt sind, so verschwindet die 
gesamte rechte Seite der Gleichung (15), und es folgt fiir alle q,: 


p"(B) + (dn — 2? cos? B) p(B) = 0. (17a) 


Da der Kern K(«, 8) auch hinsichtlich der Variablen f eine gerade Funktion 
ist, folgt ausserdem aus (9) 


p(—B) = y(+8) . (17b) 


Da in der Differentialgleichung (17a) der Parameter q nur die speziellen Werte 
_ (16d) annimmt, kénnen die nichttrivialen Lésungen mit der Eigenschaft (17b) 
_ nur die geraden Mathieuschen Funktionen erster Art ce,(6) sein. Mit q = q,(h) 
und g(a) = A, ce,(«) folgt also aus (15) bzw. (9) das Ergebnis p(6) = B,, ce, (). 
Das Verhaltnis 47 = A,/B,, ist vorlaufig noch unbekannt. Damit ist die Rich- 
tigkeit der Beziehungen (8) bewiesen unter der Voraussetzung —7 < 6 < +a. 
Das Ergebnis bleibt aber auch fiir 6 > + 2 giiltig, da das Integral in der Be- 
ziehung (8) auch fiir diese Grenzfalle existiert und den Grenzwerten stetig 
zustrebt. 
Aus (5) und (8) folgen unmittelbar die allgemeineren Beziehungen 


+7 
Cén(B) = An | Cén(x) Zo (Rk | cosa — cosB]) da (18) 
mit 
In An 
At aD Ane 


Z (x) =A Jo(x) + BN(x) und A, = 


Aus (18) folgen weiterhin die Relationen 


sin B cosa — cosB 


uaa eel. A, be | c00(e) |cosa — 00881 7 (% \cosa —cosf|) dx. (19) 
Dabei ist das Integral als Cauchyscher Hauptwert zu verstehen. 


3. Bestimmung der Eigenwerte der Integralgleichung (1) 


Um die Beziehungen (18) und (19) zur Lésung der betrachteten Integral- 
gleichungen benutzen Zu kénnen, muss man die Eigenwerte A, bzw. A, und 
4% kennen. Die Berechnung dieser Eigenwerte ist relativ schwierig und lang- 
wierig. Am unangenehmsten ist die Ermittlung der A”, und deshalb soll mit 
deren Berechnung begonnen werden. . 

An Stelle von (8) kann man mit den bereits erklarten Bezeichnungen auch 
schreiben x 

cen(B) = 2.44 | cen(e) K(x, f) de (20a) 


0 
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oder A x 


cen(B) = AN | ceq(c) R(x, B) da + 2x / ae (20b) 


0 


co 


~ 


In dem Streifen 0 < Re [a] < a gibt es nur einen Verzweigungspunkt «, = |B, 
der auf der reellen Achse liegt, da f reell sein soll (-a <[B [+ a). Von 
diesem Verzweigungspunkt fiihren wir einen Verzweigungsschnitt ber die 
reelle Achse der «Ebene nach rechts. Da die ce,(8) gerade Funktionen sind, 
darf man sich auf positive B-Werte beschranken. Unter der Voraussetzung 


0<B<a (21a) 


kann man dann die beiden in (20b) vorgesehenen Integrationswege so im 
Bogen oberhalb bzw. unterhalb des Verzweigungsschnittes von 0 nach a 
fiihren, dass tiberall auf den Wegen die Bedingung 


|cosa| > |cos B| (21b) 
erfiillt ist. Dann konvergiert wiberall auf den beiden Wegen die Entwicklung 
R(a, B) = Ny (k cosa — kcosf) = De Em Jm (k cosB) N,, (k cosa) | 1 
m=0 (21) 
Mt oe Lea eee CAO =e Be ee | 


Man kann also an Stelle von (20b) auch schreiben: 


cey(B) =A" 3" em Jpn (R cos) 


m=0 


14 % 


x / ce, (a) N,, (Rk cos a) da + ‘p Cn (x) Ny» (R cosa) a 
0 


0 
ay 


Die geraden Mathieuschen Funktionen lassen sich in folgender Form durch 
Fourier-Reihen darstellen: 


Ce, (0) me Dy 4njp COST mit:.d,.= 0 fir ++ 7=1,3,5, %..04@e 
r=0 


Damit erhalt man 


Cé,(B) == Ah Awe Jin (R cos B) Pe 
m=0 r=0 J 
4 / cosy a Ny (k cosa) da + / cosy a N,, (k cosa) da}. 
0 


0 


ca 


(23) 
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Die hier auftretenden Integrale lassen sich berechnen. Man findet mit den im 
_ Anhang kurz hergeleiteten Ergebnissen (36b): 


én(B) = 70 A SE (k cos) ) 3 da [J NO) + FeI2) NII). BA 


1= (m+r)/2  (m—r)/2 (m—1)/2 (m+r)/2 


Dabei ist nur iiber die Terme zu summieren, fiir die m+ 7 und m-+yr bzw. 
m — rv gerade sind. 

Die Doppelsummen sind fiir numerische Rechnungen unbrauchbar. Wir 
beachten deshalb, dass die Funktionswerte ce,(z/2) fiir n =0,2,4,... und 


cé,(%/2) fir nm =1,3,5,... relative Extrema sind. Zur numerischen Berech- 
nung der Eigenwerte //; eignen sich darum die folgenden Beziehungen 
ce,(2) = = 3" 2m 3) ie 14 4) n(5) fir m= 0,2,4,..., (25a) 


eal) = Broek 3) I day [Z(Z) NalZ) + F(Z) MZ | sy 
cite? =e ee | 


Ausgangspunkt fiir die Berechnung der Eigenwerte /,, ist die Beziehung 


Cop) 2 An | cen(a) Jo (Rk cosa — k cos B) dx 


Man kann nun alle Herleitungen von Beziehung (21c) bis (25b) auf diesen Fall 
iibertragen, indem man in all diesen Relationen die Neumannschen Funktio- 
nen N,, durch die Besselschen Funktionen /,, ersetzt. Damit erhalt man an 
Stelle der Ergebnisse (25a) und (25b) die folgenden: 


0e,(F eer d,, 2» [I(S)]’ fir n= 0,2,4,..., (26a) 
0e4(5) =1,2% ra Ly @, arn 15) hal) 


| (26b) 
fie) 7 3 Pee | 


Bei diesem Verfahren zur Bestimmung der Eigenwerte 4), und A, wird die 
Kenntnis der Fourier-Koeffizienten d,,, der geraden Mathieuschen Funktionen 


erster Art ce,(8) vorausgesetzt?). 


1) BE. L. Ince, Tables of the Elliptic-Cylinder Functions, Proc. Roy. Soc. Edinburgh 52, 355-423 
(1931/32). 
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Im Grenzfall k +0 gehen die ce,(8) in die Funktionen cosn B tiber. Fur; 
die zugehorigen asymptotischen Eigenwerte ergibt sich: 


1 


n 

ia ’ : Moe PRS ate. ij = Ree 
Ay = 41n (4/y k) ’ An 4 ROT re ds eo Se 

1 


> hm ) 


Ao 7 eo ee An ~0( 
mit m= n fiir n = 0, 2,4,...;m=n-+1 fir n= 1,3,5,.... 
Im Grenzfall k > 0 streben also alle #/, mit 7 = 1 mindestens wie 0(1/k%) | 
gegen Unendlich, und Aj strebt wie 1/In k gegen Null. 


4. Lésungen der beiden Integralgleichungen 


Wir kehren zuriick zu der Integralgleichung (1), wahlen aber fiir die folgen- - 
den Betrachtungen die transformierte Form (4c). Wir nehmen an, die vorge- - 
gebene Funktion G() lasse sich in eine Reihe nach Mathieuschen Funktionen | 
entwickeln. 


G(B) = 3) an cel) (0 <B <2) (7a), 


Dann folgt wegen (18) als formale Lésung der Integralgleichung (4c) der’ 
Reihenausdruck 


F(a) = 2 D" Ay ay Cen(a) (27b) 
n=0 
oder 
2 lee] 
H(a) = =— D, An Ga n(%):- (27c) 
n=(0 


Die allgemeinsten Voraussetzungen, die G() erfiillen muss, damit sowohl die 
Reihe (27a) als auch die Reihe (27b) konvergieren, sollen hier nicht untersucht 
werden. Statt dessen soll auf einen physikalisch wichtigen Grenzfall einge- 
gangen werden. 

Es wurde bereits erwahnt, dass die Integralgleichung (1) bei der Unter- 
suchung der Lichtbeugung am geraden Spalt von Bedeutung ist. Der Para- 
meter #& steht dabei in Zusammenhang mit der Spaltbreite, und zwar wird k 
um so kleiner, je schmaler der Spalt wird. Fiir den sehr schmalen Spalt erhalt 
man deshalb eine ausreichend genaue Naherungslésung, wenn man die Inte- 
gralgleichung (1) bzw. (4b) fiir den Grenzfall k > 0 list. Befindet sich ausser- 
dem die Lichtquelle in grosser Entfernung vom Spalt, so darf man 


G(B) = C = const 


setzen, und Reihe (27a) besteht im Grenzfall k +0 nur aus dem mit dem 
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Index Null gekennzeichneten Term. Deshalb braucht man in diesem Fall auch 
nur den Eigenwert Ay zu berechnen. Man erhalt aus den bereits hergeleiteten 
Grenzwerten fiir Aj und Aj und der unter (18) angegebenen Beziehung 


4 ; | 
ey ots ey ee AwhoB =. 


Dammit ergibt sich aus (27c) fiir k >0,A=1, B=7: 
UE ik aC 
2 sina In (4 i/y k) Dee lies yi — y= 21n (4i/yk)° 
Das ist genau die von SOMMERFELD!) angegebene Naherungslésung, die hier 
miihelos aus den allgemeineren Ergebnissen folgt. 
Ziir Losung der Integralgleichung (2) ist ebenso wie bei der Integralglei- 


chung (1) die Koordinatentransformation (4a) zweckmassig. Setzt man noch 
P(a) sinw = Q(«), so erhalt man dann aus (2): 


Hi) = 


|cos a — cos B| 
cosa — cos B 


W(B) = k | Qe) Z, (k\cosa —cosB|) dx. (28a) 


Wir setzen nun voraus, die vorgegebene Funktion in (28a) lasse sich wie folgt 
in eine Reihe entwickeln: 


Ley ied 
W(s) = ae on ap cenlB) - (28b) 
Dann erhalt man wegen (19) als formale Lésung von (28a) den Reihenausdruck: 
Q(x) = —2 D7 An by Cen (a) - (28c) 
nm=L 
Fs bedeutet keine Einschrankung der Allgemeinheit, dass die Reihen hier nicht 


mit dem Index Null, sondern mit dem Index Eins beginnen. Man kann namlich 
immer eine im Bereich 0 < 6 <2 konvergente Reihe folgender Art bestimmen 


DS Com Clxm(B) = 1 (29a) 
m=0 
und erhalt daraus die ebenfalls im Bereich 0 S 6 <a konvergente Reihe 
= d 
Cam —o- Cam(B) =0- (29b) 
PRL 


Ausser im Grenzfall k->0 sind alle Koeffizienten Cm in (29a) und (29b) von 
Null verschieden, da fiir & + 0 der nullte Fourier-Koeffizient dm, 9 der geraden 


1h A SOMMERFELD, Vorlesungen ber theoretische Physik, Bd. 4 (Dieterichsche Verlagsbuch- 
handlung, Wiesbaden 1950), S. 284 und 288. 
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Mathieuschen Funktionen immer von Null verschieden ist+). Daraus folgt, dass 
man im Falle k + 0 immer die Ableitung einer geraden Mathieuschen Funktion 
mit geradem Index durch die Ableitung aller tibrigen geraden Mathieuschen 
Funktionen mit geraden Indizes ausdriicken kann. Es ist zweckmassig, sich 
die Ableitung von ce,(8) durch die Ableitung der iibrigen ce,,,(8) ausgedriickt 
zu denken, um bei der Behandlung des Grenzfalles k = 0 keine Schwierigkeit 
zu haben. Deshalb fehlen in den Reihen (28b) und (28c) die Terme zum 
Index Null. 

Die Lésung (28c) der Integralgleichung (28a) ist im allgemeinen nicht 
vollstindig. Wenn namlich die Eigenwerte A, so beschaffen sind, dass die mit 
den Koeffizienten c, der Reihe (29a) gebildete Reihe 


dD, 4am Cam lem(%) = Qo(x) im Bereich 05a (29¢) 
m=0 


konvergiert, so stellt Qo(a) wegen (29b) und (19) eine Lésung der homogenen 
Integralgleichung (28a), also eine Lésung fiir den Sonderfall W(f) =0 dar. 
Die allgemeine Lésung ist dann also 


Ola) = =2 x An bn Ce (@t) + © Qole) (30) 


wobei der konstante Faktor C beliebig gewahlt werden darf. 

_ Auf die Ermittlung der Voraussetzungen, unter denen Reihe (29c) konver- 
giert, soll in dieser Arbeit nicht eingegangen werden. Als notwendige Bedingung 
ergibt sich, dass der in der Kernfunktion 


ZAM) a fi AB NARS 


auftretende Faktor B von Null verschieden sein muss. Hinreichende Bedingun- 
gen fiir die Existenz der Reihen (28b) und (29c) gibt SCHUBERT?) an. Wenn 
B von Null verschieden ist, diirften hier die Verhaltnisse ahnlich liegen wie in 
dem eingehend untersuchten Sonderfall k > 0. 

Zum Abschluss wollen wir einige der mitgeteilten Ergebnisse an dem ein- 
fachen Sonderfall & - 0 beleuchten. Damit in diesem Grenzfall das Integral 
der Integralgleichung nicht identisch Null wird, muss der in der Kernfunktion 
auftretende Faktor B von Null verschieden sein. Im Grenzfall k +0 nimmt 
(28a) die Form an 


7 cosa — cos B 


wis ees / oie eras 
0 


1) E. L. Incz, a.a.O., S. 358. 
2) H. Scuusert, Uber eine lineare Integrodifferentialgleichung mit Zusatzkern, Ber. Verh. Sachs. 
Akad. Wiss. Leipzig, Math.-naturw. Kl., 97, Heft 7 (1950). 
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Die geraden Mathieuschen Funktionen ce,(f) sind hier zu ersetzen durch 
die Funktionen cos 6 mit n = 0, 1, 2, .... Reihe (29a) entartet zu einem einzi- 
gen Term mit dem Index 0, und man findet cy = 1. Reihe (29b) entartet zu 
der trivialen Aussage (d/d8) cosn B= 0 fiir »=0. Dementsprechend ist 
Qo(«) = My. Mit den bereits angegebenen Grenzwerten fiir 2” nimmt also die 
Lésung.(30b) die Form an 


1 oo 
Oa) = BE DS nb, cosn «% + const , 


n=1 
wenn sich die vorgegebene Funktion in eine konvergente Reihe der Form 


1 «< : 
W(B) = = sin B Ee nN by, sin 7 p 


entwickeln lasst. Setzt man 2 B = —1 und 1b, =d,, so hat man damit die 
bekannten Ergebnisse aus der Tragfliigeltheorie, die sich bei uns als Grenzfall 
von allgemeineren Beziehungen ergeben. 


5. Anhang: Berechnung einiger Integrale 


Ausgangspunkt der Berechnung ist die bekannte Beziehung 


7/2 


i JA q COS %) COSY & da = > Jere 2) Jus (uu es =1) (31) 
0 


Wir wollen nicht von 0 bis 2/2, sondern von 0 bis z integrieren. Da aber fiir 
nicht ganze y-Werte der Punkt « = 7/2 Verzweigungspunkt ist, miissen wir 
fiir einen solchen Integrationsweg angeben, wie der Verzweigungspunkt um- 
gangen werden soll. Wir definieren zu dem Zweck zwei Integrale. Durch das 
Zeichen - bzw. v werde ein Integrationsweg beschrieben, der, von 0 ausgehend 
oberhalb bzw. unterhalb am Verzweigungspunkt 2/2 vorbei nach z geftihrt 
wird. Es sei also 


ID) py = [ JAZ qgcosa) cosnadx mit n= OVI 9S Press (32a) 
0 


fy 
Aus den Umlaufrelationen fiir Bessel-Funktionen und der Periodizitat von 
cos n « folgt unter der Voraussetzung uw = —1: 
7/2 


D,, n=l[l+(-1"" efeny | JAZ q cosa) COs n % dx 
0 
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oder wegen (31) 


Din => + (1 el] Jum @ Juma - (32b) 


Der Integralausdruck (32a) existiert fiir beliebige reelle ~-Werte und stellt eine 
analytische Funktion von y dar, sofern g + 0 ist. Deshalb gilt das Ergebnis 
(32b) nicht nur unter der Voraussetzung « = —1, sondern fiir beliebige reetle 
Werte von wu. 

Fiir das zweite Integral 


D,, = | J,A2qcosa) cosnadx mit n=0,1,2,3,... (33a) 
v 
findet man auf Grund derselben Uberlegung das fiir beliebige reelle ~-Werte 
giltige Ergebnis 
Dan = ZL + 1)" 698] Jounal) Jin—nyold) - (33b) 


Durch Addition von (32) und (33) folgt 
Ds + De = [ F,A2 g cosa) cosm « dx + / .. ee | 
é a (34) 
= a{1+ (=1)? COS 7] Tins ny 2(q pe oe ny2(9) | 


fiir n = 0,1, 2, 3,... und beliebig reelle u-Werte. 
Im Anschluss an (32) und (33) definieren wir nun noch die folgenden 
beiden Integrale: 


Ee = f Nal2q cose cos nada, Em» =f de (35) 
0 0 


mui 0,1,2,3;... und ee. 0,122. 5.5.0. 
Auf Grund der bekannten Beziehung 


Bae | eu Ose 
oe , By . a on Yrs It m 


kann man die Integrale (35) durch Differentiation aus den Integralen (32) und 
(33) herleiten. Man findet: 
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Daraus erhalt man schliesslich nach einer Reihe von Umformungen 
oe am m 
aa a) ( 18 ah J m+n) 29) Jim — nyo) a oe (—1) J(-m+ny2(Q) J(~m—ny2(Q) 


ott > [1 ag (eke [Tim +-my2(2) Non —ny2(2) 75 Jot Leno) Nom +ny 22) | 
(36a) 


Hieraus erhalt man die Beziehung fiir es mn» ndem man iiberall + durch —7z 
ersetzt. Die Addition der beiden Ausdriicke ergibt: 


& 
= 
= 
a 
gy 
I 


oe HENS (2 q cosa) cosmada + |... de 
0 0 


I 


= [1 = ig auees] [Jin +-nyj2(Z) Non —ny2(9) aF Jon —ny2) Nin +ny2(9)]- 
(36 b) 


Damit haben wir alle in den vorhergehenden Abschnitten benutzten Integrale 
hergeleitet. 


Summary 
The author investigates the integral equation of the first kind with the kernel 
A Jo[k (cosa — cos B)] + BNj(k |cos« — cosB)) . 


This kernel is shown to have the even Mathieu functions as eigenfunctions. A 
method suitable for numerical computation of the corresponding eigenvalues is 
given, and the solution of the integral equation from the obtained results is 
outlined. By partial differentiation of the kernel with respect to B an equation 
is obtained which may be solved using the previous results. Both equations 
appear either in the theory of diffraction of sound or light by slit of infinite 
length or in the theory of a thin airfoil oscillating in a compressible fluid. 


Résumé 


On présente une étude de l’équation intégrale de premiére espéce ayant pour 
u 
ie A Jolk (cos« — cos B)] + BNg(k |cosa — cos B)) . 
On montre que ce noyau a pour fonctions propres les fonctions de MATHIEU paires. 
On donne une méthode de calcul des valeurs propres associées appropriée a un 
calcul numérique, et on indique le moyen de parvenir ala solution de l’équation 
intégrale en utilisant les résultats obtenus. Par differentiation partielle du noyau 
par rapport a B on obtient une équation apparentee a l’équation considérée et 
dont la solution se déduit d’une maniére relativement simple de la solution de 
cette derniére. Les deux équations interviennent dans le traitement de la diffrac- 
tion des ondes sonores ou lumineuses par une fente infiniment longue et dans le 
calcul des forces aérodynamiques mises en jeu par un profil mince oscillant har- 
moniquement dans un milieu compressible. 


(Eingegangen: 20. Mai 1952.) 
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Sur la thermodynamique des processus irréversibles 


(Second mémoire) 


Par KyrIL_E Poporr, Sofia?) 


Introduction 


Dans un premier mémoire?) nous nous sommes occupés de la théorie des 
processus irréversibles telle qu’elle a été esquiss¢ée par ONSAGER®), développée 
par Casimir4), DE Groot’) et d’autres et fondée sur la symétrie du tenseur L 
des relations phénoménologiques 


if —— Qe i xX, ; avec Lix = Ly; ? (1) 
k 


ou les J; représentent des flux et les X, des forces, c’est-a-dire les dérivées par 
rapport au temps ¢ des x; et les dérivées partielles par rapport 4 x; de —AS 
respectivement, AS étant l’expression de l’entropie, autour de l’état stable 
d’un systéme adiabatiquement isolé, réduit a 


2 1 
AS = = oD Bea Xp, avec Sir = Sxi- (2) 


Dans notre premier mémoire nous avons montré que dans ce cas les relations 
phénoménologiques représentent des intégrales premiéres du systéme d’équa- 
tions différentielles 


d?x; ~ 

m4 , a 6 

aa = avec Ay = —- ee (3) 
que la symétrie du tenseur L;;, est une conséquence logique de (3), dans ce cas 
réduit, et que les L,;, sont des constantes. 


Dans l'étude qui va suivre nous considérons AS sous sa forme la plus 
générale et nous ¢tablissons que les relations phénoménologiques découlent des 


1) Institut de Mathématiques 4 l'Université de Sofia. 

*) Kyritte Poporr, Sur la thermodynamique des processus irréversibles, ZAMP 3, 42-51 (1951). 

3) L. Onsacer, Reciprocal Relations in Irreversible Processes, Phys. Rev. 37, 405-426 (1931)5 
38, 2265-2279 (1931). 

4) H. B.C. Casimir, On Onsager’s Principle of Microscopic Reversibility, Rev. mod. Phys. 17, 
348-350 (1945); Philips Res. Rep. 7, 185 (1946). 

5) S. R. pe Groot, Thermodynamics of Irreversible Processes (North-Holland Publishing Com- 
pany, Amsterdam 1951), 

8) Comme il a été mis dans les Errata du Fascicule 3, page 237, de ZAMP, les calculs dans notre 
premier mémoire sont faits avec les équations (3) ci-dessus. Au lieu de AS lire —AS dans les for- 
mules (4), (5), (10), (11) du premier mémoire. 
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équations (3), mais dans ce cas la symétrie du tenseur L;;,ne se manifeste que 

_dans les termes de la premiére approximation. Pour arriver a ces résultats, il 
nous a fallu avoir recours a la théorie analytique des équations différentielles 
établie par HENRI Poincaré!) et développée par EMILE PICARD?). Aprés avoir 

_démontré que les équations (3) admettent des intégrales s’annulant pour ¢ = +co 
et établi le caractére analytique de ces intégrales, nous en déduisons l’ordre de 
grandeur de la différence L;, — L,, dans le cas général. 


Pose du probléme et étude des équations différentielles correspondantes 


Soient &), &,..., €,, les paramétres définissant l'état du systéme, s(&,, &,..., €,) 

Ventropie du systeme a2 un moment donné, &9, &,..., €9 les valeurs de ces 
parametres, correspondant a l’état stable, défini par le maximum de S. On a, 
dans le cas de deux paramétres &,, &, en posant 


Be Ey Gy yn Vie, 
; 1 (02S as as 4 
BS 5 (E77 62) — S(T, 62) = 5 rs xe 2 ort oye ) (4) 
1 /0%S DEST oe ) 
| 3! (See RG Ox? ie UE Ox dy? ~ y “Oy Ae Y 


les dérivées partielles se rapportant aux valeurs x = y = 0. 
Ainsi, dans le cas de deux variables, les équations différentielles (3) peuvent 
étre écrites sous la forme 


ge Xa aet bytes 2payteyt-, | 

: (5) 
2 2 
Cee ae Loy+fxu*+2exythy?+te--, | 


ota, b, c, ... remplacent des dérivées partielles de S, pour x = y = 0, multipliées 
par des constantes convenables, a «* + 2 bx y+ cy? étant une forme quadra- 
tique définie positivement et oti, par conséquent, l’on a 


ped 6. ae SOG BOE (6) 
Le systéme (5) est équivalant au systeme ) 
et a l’équation aux dérivées partielles 
oF gy Ex 4+S H+ 5 ¥=0 (3) 


1) HENRI Poincare Sur les propriétés des fonctions définies par les équations aux différences pur- 
2 / 


tielles, These de doctorat (Paris 1879). We 
2) EMILE PICARD, Traité d’analyse, tome 8 (Gauthier-Villars, Paris 1908). 
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Pour étudier le caractére analytique des intégrales x, €, y, 1 du systéme (7) 
s’annulant simultanément, ou de l’équation équivalante (8), faisons le change- 
ment de variables défini par 


U = Oy E+ Oy X + yg N+ Mg Y, W= Ay S + Ogg X + Ogg) + Age Y | 


(9) 
U = Og E + Ogg X + Og M+ Oleg Y, = Og S + Oye X + Ogg + Maa Y - 
On a ainsi 
OF OF OF OF OF 
OE Ou AE Gy MET Gy Mat Tog Mat? 
OF OF OF OF OF 
Oe Gu MT Gy MET Oy Mee Foz Mae ia 
OF OF OF OF OF ( ) 
On Ou “8 T Gy MT Gy Mas Tog Mea? 
OF OF OF OF OF 
Oy) bu Pe Cad ae Cea Ona 
En portant ces valeurs des dérivées de F dans l’équation (8) on obtient 
OF es 
Wy (041 X -+ aye E+ a3 ¥ + O44 H) 
OF , =a 
hase (Gq) AX + O99 5 + 93 Y + ogy H) 
(11), 


OF 


sali (ag, X + ogo SF + %g3 Y + X54 H) 


OF ; = 
Te (Oq1 X + Og. F + O43 Y + O14, H) = 0. 


Choisissons les «,;, de fagon qu’on ait, aprés la réduction, 
OF OF OF OF 
me UE ES QV) +Oe (wt) tS (at) =0, (12) 


U 


les points... remplagant des termes de second degré et de degrés supérieurs 
a2en wu, v, w, 2. A l’équation (12) correspondra ainsi le systéme d’équations | 
différentielles 
au e dv = dw dz me a (13) 
V1 U -t- +++ 


YoU +o ¥5W + TE ae 

Pour cela il faut avoir 
Oy Ay + Oye B+ Og Vi + yg HH =7y (Gy E+ Oy % + 497 + 4), 
gy Xy + Ogg E+ dys Vy + Ong H = 1q 0 = Me (tg, & + Ogg % + O93 9 + O04) 5 (14) 
gy Xy + Ogg E+ Otg5 Yi + Ogg H = 15 W = 1p (X51 § + Og ¥ + Orgy 9 + Oya ¥) , 
gy Xy + Ogg SF + Ogg Vy + otgg H = ry 2 = Hq (tg + Orgy % + O43) + Hayy) , 


ou X, et Y, désignent les X et Y réduits a leurs termes du premier degré. 
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On tire d’ici, pour déterminer les 7, le systéme 
hip @ — hee %e +30 +a4°-0=0, | 
— O97, Php 01 + oe? OF a4 0=0, 
(15) 
0 + hig: O+ aj3¢ — bir, =O, 
31° O+ aj: O — a,57; ain Chen © L=@, 


Pour que ces Equations, homogenes en «; ,, Solent compatibles, il faut qu’on ait 


GD 0 | 
—y, 1 0 0 | 
ay (16) 
b 0 Co er, 
0) OL = | 
ce qui donne pour déterminer 7; l’équation 
yi— rv? (atc) t+ac—b*=0 (16 bis) 


Comme nous l’avons vu, dans notre premier mémoire, ces racines sont 
réelles et l'on a 


ae Oi 7 ee (17) 


Avec ces valeurs de 7;, les équations (15) donnent les «;;,. En posant «,; = 1, 
‘lon obtient ainsi 


v2 — a ar 
CCS Ou ge Bi Gc (tae a Gs ae, bg = 17% — Yi» (18) 


a;, Bir 7; étant les valeurs de ces quantités d’aprés notre premier mémoire cite. 
Les intégrales des équations (13) tendant vers zéro pour des valeurs positives 
de ¢ sont, comme on le sait des travaux de PoIncarE et de PicarD cités, des 


fonctions holomorphes de 
IRS HOS, 


K, et K, étant des constantes arbitraire assez petites. Elles peuvent étre déve- 
loppées en séries suivant les puissances entiéres et positives de Ky eaieiere. 
convergentes dans l’intervalle (0, + oo) de ¢t et tendant vers zéro pour ? = + 00. 
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Revenons aux équations (7) et au changement de variables (9). Le déter- 
minant des coefficients de x, vy, &, 7 des relations (9) étant 
Le Ss a | 
1 / 
. | = 41 14 (ty — %4)? + 0, 
1 fs % Ys / 
1 


ta hae a 


ces équations déterminent x, y, &,7 comme fonctions linéaires des variables 
u,v, W, Z, s'annulant avec eux. Par conséquent les intégrales x, y, €,7 du 
systéme (7), s’annulant pour ¢ = + oo, peuvent étre développées de méme en 
séries procédant suivant les puissances entiéres et positives de K, e’*’, Kye”, 
convergentes pour ¢ = 0. On a ainsi 
Yes &, ies erst -}. Gs a ent ae G. RK ez! A Gs he K, eltztr)t a - Kie 2Qr,t “oa 
y = G, K, &' + G, K, e' + G, K3 ce?’ + G, Ka K,e*"* + G, K? ee free 
F (19) 
d’ou l’on tire 
d*x 


Saw =O Kye + C, Kyrie + C, Kp 4g ot! 
|. G, Ke KG (ro afte 74)% est r,t ae CG Kk 4 ‘A e2nt Bar 


et une expression analogue pour d*y/dé*, ot les coefficients C; sont remplacés — 
par G;. 

Pour trouver les valeurs de ces coefficients, il n’y a que former les expres- — 
sions de X et de Y avec ces valeurs de x et y et employer la méthode des 
coefficients indéterminés. On obtient ainsi: 


X=ax+by+ex*+2fxy+gy?t+--: 
= K, (aC, +b Gi) &! + Ky (aC, + 6G) e* 

+ K3 (aC, +0 G+ eC? + 27C, G, + g G3) et" 

+ Ky K,(aC,y+6G,+2eC,C,+ 27 (C,G,+C,G,) +22G, G,) e' 

+ KG (aCy+bG,+eC3+27C, G+ g G3) er +..., (20) 
Y = Ky (bC, +c G,) e* + K, (bC,+¢G,) & 

+ K}(bCs+¢G,+ $C? +2gC,G, +h Ge" 

+ K, K,[6C,+cG,+2/C,C,+2¢(C,G,+C,G Dende | 

+ KG (C3+cG5+ fC}+2gC,G,+h G3) e+ ..., 


La comparaison des coefficients du méme degré en e” et e” dans les ex- 
pressions de X et de d?x/dt® ainsi que dans celles de Y et d?y/dt? conduit tout 
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d’abord aux équations suivantes: 


Ci (78 —a)—bG,=0, b C, + (72 — c) G, =0; | 
(21) 
Crp) = 02 OG, Se) G, = 0.5 | 
Les deux premiéres de ces équations donnent pour G, 
v3— a bC 
G, — Ce ia et Gy = Seer ) (22) 


ces deux valeurs de G, sont égales d’aprés la relation 
(3 — a) (8 — 0) = BY, 
servant a déterminer la racine 7,. Comme d’autre part 


v3 — a 
Tel = A , 


on a G,= Cy, . De méme les deux valeurs de G, tirées du second systéme 
sont égales, et on a G, = Cy a. 
La comparaison des coefficients de e?*, e?”", etn)? |. conduit aux 
systémes d’équations suivants pour la détermination de C;; G;: 
Ca (472 —a) —G,b=eC?4+27C,6,+ 2G, 
0,04 G,(44—¢c)={CG+2¢C,G, +4 Gi; 
Cl Gerrs)? — 2] G,6=2eC,C,+27 (C,G,+ CG) +2¢G, &, 
ea. Dt 08 4 74)" c]G,=27C, C,+ 2 ¢ (Cy G+ Ce G) + 2h G, Gy; 
Ce (472 — a) G,bs¢C24+2/C,G,+ 2G, 
= C0 G, (474 — 0) = }/C2+2¢C,G,+h G3; etc. 


Il n’y a pas lieu de s’occuper ici de la convergence des séries ainsi obtenues, 
la théorie générale de ces équations nous assurant d’avance cette convergence. 


Dans le cas général les relations phénoménologiques découlent des 
équations (5) 
Pour montrer que, dans le cas général, les relations phénoménologiques 
sont des conséquences logiques de 


1 a(4s) d¥y _ y _ _ (AS) (5 bis) 


ro ee ea ae oy? 
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mettons les intégrales (19), en ne gardant que des termes jusqu’au second 
ordre en K,, K,, sous la forme 


, Y 1 Ui eg 
n=O" Ky C, (L+ Ch Kye! + 5G Kye) 


Pea eyes (1 Oe Rie ae 


, 2 1 , ay 
y =e Ky G, (1+ Gy Kye + > G Kye) 


+ eK, G, (14+ > G, Kye + Gs Ky é), 


ou l’on a mis 


On a méme 
x = eK, Cr ( tick 2 Cy Ka ten ake, e) Hin et Ky CK, tee 
4+! Ky Cor, (1 + > C, Ke +20, K, zt) oo" EC, KC ame 
y’ =e" Ky Gyr (1 + 2G, Kad 4+ 2 GK, o) +5 et KG, KG re" ; 


+ eK. Got, (1 of pes K, &' + 2G, Ky ) 4. = e' KG, KG, t, &™. 
(25) } 
On a d’autre part G, = Cy %, Gy = Cy «4. Avec cela nous pouvons écrire 


xm et KOC (le A ee A Cyl ea), | (26) 
y =e Ky Cy ag (1+ 0) te Ky Cyay (1 +2), | 
x= oY KC (Ld) POOR, Cyl +a), one 


y =O Ky Cyryeg (lta) te Ki Cam (Ltm), | 


ott A, M4, 0, %, Ay, fy, Q1,%, Sont de premier ordre par rapport a K, et K,. 
En substituant a e’, e’' dans les derniéres équations leurs valeurs tirées 
des deux premicres, l’on obtient 


x6 = [ry oq (1 +2) (L+A,) — 14% (1+ @) (1+ )] 
+ y [ra (+A) (1+ ay) — 72 (1+ w) (1 +4,)], 

YO= [Fo Hp 4 (I +7) (1+ @1) — 74 Hy & (1 + 0) (1+ 7,)] 
+ y [rq aq (1 + A) (L +g) — 72% (1+ w) (1+ @,)], 


(28) 
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ot. lon a 
b= (144) Ata) y—G+p) (+o) m.- (29) 


Pour obtenir les relations phénoménologiques il n’y a qu’a exprimer x et y 
au moyen de X et Y. On a pour cela 


Aaa leave SY = oe IT 2 4) IN (30) 
avec 


L=a+ex+fy, M=b+fut+egy, N=c+exthy. 
D’ici l’on tire 
xA=XN—-YM, yA=—-XM+YL (31) 


avec 


A=LN—M?. 


Avec ces valeurs de x et y, les relations (28) deviennent 


Adz’ = X{N[reoq (1+) (L +A) — 74% (1 + @) (1+ My) | 
—M[n (i +4) +m) — (+a) + a)] 
4 V(L[n (tA +m) —n ty) +a) 
— M [ro oq (1 +2) (1 + Ax) — 74% (1 + @) (1+ u)]} 
=AX+BY 
Ady’ = XN [ra og 4 (1 +2) (1+ 01) — 74% 0a (1 + @ )(L + ™))] 
—M [rg (14 A) (1 + ty) — 1% % (1 + p) (1 + 01) |} 
+ V{L [re og (1 +A) (1 + ata) = 72% (1+ w) (1 + 4) 
— M [1re %» %4 (1 + 2) (1 + 01) — 14 % % (1 + @) (1 + 7,)]} 
=DX+CY. | 


(32) 


Ce sont les relations phénoménologiques. Nous allons montrer que dans le 
cas général que nous considérons ici, la tenseur Ej, est presque ee 
c’est-A-dire que les coefficients B et D ne différent que par des termes qui con- 
tiennent K, ou K, en facteur. En effet, en ne gardant dans les expressions de 


ces coefficients que des termes ne contenant pas K, et K, d'une maniére directe 


ou indirecte, il faut mettre 


aa Mah Neat 
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et par conséquent 
b 5 2 
B= a (14— 1%) — B (1 %q — 14%) = 4 (14 — 72) — ae (i= a7, 05 a)| (33) 
= 1214 (Yo — 74) 


et 
D = C (Yq % 4 — 14 4g %q) — B (Fy Hy — 7y Mp) - (34) 


Mais ¢ a = & 7% — b, Cu = Hy 77 — 0D et avec ces valeurs de ca; ona 
D = Oy %4 (74 — 1) 1214 - 


D’autre part 


(3—a) (5-4) Ar -—alre+ r+ a 1 
Xs Hy -s se? b2 ae — = b2 =. = 
et par conséquent 


Ce quimontre qu’on ait Lj, = L,, aux termes du premier ordre en A, et Ky pres. 


Zusammenfassung 


In einer ersten Abhandlung haben wir uns mit der Theorie der irreversiblen : 


Prozesse beschaftigt, die von ONSAGER und CASIMIR entworfen ist und welche 
auf die Symmetrie des Tensors L;, in der phanomenologischen Relation 


vo 2 Lin Xx 


beruht. Hier bedeuten J; die «Flux» und X, die « Krafte», das heisst die Ableitung 
von #; in bezug auf die Zeit ¢ bzw. die partielle Ableitung von —AS in bezug 
auf ¥,, wo AS den reduzierten Ausdruck 


1 : 
AS = ~ a & bik % Ve, Mit Sin = Sri» 


der Entropie in der Nahe des Stabilitétszustands eines adiabatisch isolierten 
Systems bedeutet. 

Um die Symmetrie des Tensors L;;, nachzuweisen, gehen ONSAGER und CasI- 
MIR von der Theorie der Wahrscheinlichkeitsrechnung aus. In unserer ersten Ab- 
handlung haben wir gezeigt, dass die Symmetrie des Tensors L;, eine Folge des 
Gleichungssystems 

d*x; xX : 

aa = Xi Hae Me Pa a 
ist und dass die L,;; konstant sind. In der vorliegenden Abhandlung zeigen wir, 
indem wir von der analytischen Theorie der Differentialgleichungen ausgehen 
und die Glieder von drittem und héherem Grade von 4S beriicksichtigen, dass 
in diesem Falle die Hauptteile der Ausdriicke von L;;, symmetrisch sind. 


(Regu le 15 avril 1952.) 
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Rechengerate mit linearen Potentiometern’) 


Von AMBROS P. SPEISER, Ziirich?) 


Einleitung 


Es ist bekannt, dass sich Netzwerke mit variablen ohmschen Widerstanden 
unter gewissen Voraussetzungen als Rechengerate verwenden lassen, welche 
sich durch besondere Einfachheit auszeichnen. Unter Verwendung von linearen 
Potentiometern kénnen Gerdte zur Bildung von skalaren Produkten (und 
damit zur Lésung gewisser linearer Gleichungssysteme durch Iteration) gebaut 
werden [1]8). Mit Stufenwiderstanden, die durch entsprechend vielschrittige 
Stufenschalter betatigt werden, sind bemerkenswerte Rechengerate zur Aus- 
wertung komplizierter trigonometrischer Ausdriicke entwickelt worden [3], [4], 
[5]. Potentiometer kénnen mit einer nichtlinearen Widerstandskurve gewickelt 
und so zur Erzeugung von Funktionen eingesetzt werden [6], [7]. - Dagegen 
kénnen Differentiation und Integration mit ohmschen Netzwerken nicht aus- 
gefiihrt werden. 

Die vorliegende Arbeit zeigt, wie Rechengerate konstruiert werden k6nnen, 
deren Funktion iiber die blosse Bildung von skalaren Produkten hinausgeht, 
in denen aber trotzdem neben ohmschen Widerstanden nur lineare Potentio- 


meter verwendet sind. 


1. Die elementaren Rechenoperationen 


In diesem Abschnitt werden die hier zu verwendenden Schaltungen fiir die 
elementaren Rechenoperationen rekapituliert [1], [6], [7]. Die zu verarbeiten- 
den Variabeln sind als Drehwinkel, Spannungen oder Strome gegeben. 

Figur 1 zeigt eine Schaltung zur M ultiplikation der Grossen x und y. Hier 
(wie auch im folgenden Abschnitt) kennzeichnen x und y die Stellung der 
Biirste auf einem Potentiometer. Diese Variabien legen zwischen 0 und 1, 
wobei der Wert 0 der geerdeten Seite des Potentiometers zugeordnet ist. 
U ist die Spannung der speisenden Quelle; alle Spannungen sind gegen Erde 
gemessen. 

Es ist ersichtlich, dass nur solche Variable miteinander multipliziert werden 
k6nnen, die als Drehwinkel gegeben sind. Das Resultat erscheint als Spannung 


1) Eine erweiterte Fassung dieser Arbeit befindet sich an der Bibliothek der ETH. in Ziirich [8]. 
Dort sind weitere Spezialschaltungen gegeben sowie eine Erlauterung der nomographischen lurven- 
anpassung; ferner sind vollstandige Rechengerate beschrieben. 

2) Institut fiir angewandte Mathematik an der ETH. ' a - pie a” 

3) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis aul Seite 459. 
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und ist nur richtig, wenn die Konstante « gross gegeniiber 1 ist; infolge der 
Belastung des ersten Potentiometers durch das zweite weicht die erhaltene 
Spannung « um Au vom richtigen Wert ab. Die grésste relative Abweichung 
entsteht bei « = 2/3 und hat fiir y = 1 den Wert 

Au 0,15 


re ao 


U 
x 
Rk 
eR SG lig ie xy 
Fig. 1 Fig. 2 
Multiplikation. Addition. 


(Diese Beziehung gilt fiir « > 1.) Mit « = 20 wird also die relative Genauigkeit 
besser als 19%. — Es besteht aber auch die Méglichkeit, die Skala fiir % so ein- 
zuteilen, dass der durch die Belastung entstehende Fehler ausgeglichen wird. 
Dadurch fallt die Bedingung « >> 1 dahin. 

Die Addition von Spannungen kann so geschehen, dass diese in Serie ge- 
schaltet werden. Dazu sind aber unabhangige Stromquellen nétig (etwa meh- 
rere Batterien oder mehrere Sekundarwicklungen auf einem Transformator). 
Da hier Gerate betrachtet werden, die von einer einzigen Batterie gespeist sein 
sollen, fallt dieses Verfahren ausser Betracht, und es wird die Schaltung von 
ligur 2 verwendet. Um die Summe wu, + u, zu bilden, wird die Biirste z so 
lange verstellt, bis der Strom J, der durch das Messinstrument fliesst, ver- 
schwindet. Dann liegt die mit S.S. (Sammelschiene) bezeichnete Leitung auf 
Erdpotential, und es gilt: 

9 VAG] 
Sisk mart ite (1) 


oder, falls die Widerstande R,, R, und Ro, die als Summationswiderstinde be- 
zeichnet werden, alle gleich sind: 


Us 
RR, a R. 


zU = —(u,+ u,) . 
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(Ahnlich wie bei der Multiplikation gilt auch hier die Einschrankung R, > R.) 
Die gesuchte Summe kann also an der Skala x abgelesen werden. Dies ist eine 
«implizite Losung», welche durch einen als Nullabgleich bezeichneten Prozess 
gewonnen wird; die Gleichung (1) kann ohne Anderung der Schaltung auch 
nach u, oder u, aufgelost werden, falls die anderen Gréssen gegeben sind. 

Die gezeigte Anordnung ist nicht auf zwei Summanden beschrankt ; vielmehr 
konnen auf die Sammelschiene beliebig viele Summationswiderstande geschal- 
tet werden, wobei der Nullabgleich natiirlich nur einmal auszufiihren ist. Sind 
die Spannungen ,, w#, ... die Resultate von Multiplikationen gemass Figur 1, 
so ist das Resultat ein skalares Produkt. (Dabei darf jedoch die Belastung des 
zweiten Potentiometers in Figur1 nicht ausser acht gelassen werden.) Durch 
eine entsprechende Schaltung lasst sich auch x y = z bilden; da es sich auch 
hier um ein implizites Lésungsverfahren handelt, kann diese Beziehung als 
x = —z/y geschrieben werden, somit wird die Operation des Dividierens aus- 
gefiihrt. 

Das Instrument zur Messung von / muss so empfindlich sein, dass der Null- 
abgleich mit hinreichender Genauigkeit erfolgen kann. Um aber Uberlastungen 
bei den zuerst fliessenden, hohen Fehlerstro6men zu vermeiden, soll die volle 
Empfindlichkeit erst eingeschaltet werden, wenn der Abgleich annahernd 
richtig ausgefiihrt ist. An Stelle des Instrumentes kann nattirlich auch eine 
Servosteuerung mit Motor treten, welche den Nullabgleich automatisch aus- 
fiihrt. Dadurch geht aber die Einfachheit und Anspruchslosigkeit des Gerates 
verloren. 


2. Die Erzeugung von Funktionen einer Variablen 


Die einfachste Art der Erzeugung einer Funktion ist eine entsprechende 
Einteilung der Skala, an der die unabhangige Variable eingestellt wird. Dieses 
Verfahren versagt jedoch, sobald von derselben Variablen mehrere verschie- 
dene Funktionen gebildet werden sollen. In solchen Fallen muss von der Fahig- 
keit linearer Potentiometer zur Bildung von Funktionen Gebrauch gemacht 
werden. Es gilt folgender Satz: 

Gegeben eine Schaltung (Figur 3), die von der konstanten Spannung U ge- 
speist ist und die auf einer gemeinsamen Achse mit dem Drehwinkel x (0 S x £1) 
endlich viele lineare Potentiometer, ferner feste Widerstdnde enthdlt. Dann ist die 
Spannung u(x) an einem beliebigen Punkt dieser Schaltung eine rationale Funk- 
tion von x. Diese Funktion hat im gegebenen Bereich keine Pole, und Nullstellen 
hinnen nur am Rande des Bereiches vorkommen. 

Dieser Satz geht aus den Gesetzen fiir die Serie- und Parallelschaltung von 
Widerstanden hervor. — Die Mannigfaltigkeit der méglichen Schaltungen ist 
dabei ausserordentlich gross. Beschrankt man aber die Zahl der Potentiometer 
auf zwei, so lassen sich die Schaltungen klassifizieren [11]. Es bleiben dann 
gewisse Typen von rationalen Funktionen, die elektrisch erzeugt werden, und 
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mit diesen lassen sich viele in der Praxis vorgegebene Funktionen — seien diese 
nun analytisch oder als Kurven gegeben — approximieren. Es wird hier vom 
wichtigen Prinzip Gebrauch gemacht, dass eine verlangte Funktion oder Ope- 
ration durch eine andere, prinzipiell davon abweichende ersetzt wird, so dass 


U 
x 
= u(x) 
U ufx) R 
ak 
Fig. 3 Fig. 4 
Schaltung zur Erzeugung von Funktionen. Das belastete Potentiometer. 


gd 


Fig. 5 


Funktionen des belasteten Potentiometers. 


eine exakte Wiedergabe selbst mit fehlerfreien technischen Mitteln nicht mehr 
moglich ist. Da aber die verwendeten Potentiometer mit Ungenauigkeiten be- 
haftet sind, geniigen approximative Darstellungsarten fiir vorgegebene Funk- 
tionen, sofern deren Abweichung in der gleichen Gréssenordnung wie die Fehler 
der Potentiometer liegen. Es ist zu bemerken, dass unter Erhdhung des schal- 
tungstechnischen Aufwandes die Abweichungen in dieser Darstellung beliebig 
klein gemacht werden kénnen. 
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Besonders ausgiebig untersucht wurde die Schaltung des belasteten Poten- 
tiometers, die in Figur 4 gezeigt ist. Die Spannung an der Ausgangsklemme ist 
folgende Funktion von «x: 


ey 


aN 
(UE g(x) ee Lp 


Der Verlauf von g(x) ist fiir verschiedene « in Figur 5 gezeigt. Zur Darstellung 
einer vorgegebenen Funktion kann eine Kurve aus dieser Schar ausgewahlt 


+Y 


uf(g)— 
U sin p 


Fig. 6 


Schaltung zur Darstellung der Funktion sing. 


werden. Es kann auch aus dem gezeichneten Bereich nur ein Ausschnitt ver- 
wendet werden; zu diesem Zweck ist in Serie zum Potentiometer in Figur 6 
oben und unten je ein fester Widerstand zu schalten. Diese Widerstande mogen 
die Werte B R und y R haben. Indem nun «, 6 und y frei wahlbar sind, haben 
wir zur Anpassung eine dreiparametrige Schar von Kurven zur Verfiigung. 
Aus Figur 5 ist ersichtlich, dass man damit eine grosse Mannigfaltigkeit von yer 
schiedenen Formen zur Verfiigung hat. Die Kurven sind nach unten gekriimmt ; 
durch Bildung der Umkehrfunktion entsteht eine entsprechende, nach oben 
il char. Jee 
Sc tcchwiovigeet ist die, fiir eine vorgegebene Kurve dasjemige Tripel 
der Parameter zu finden, das die beste Annaherung gewahrleistet. Zahlreiche 
Versuche haben zu folgenden Erkenntnissen gefiihrt : Pare Naherungslosungen 
(bis zu einer Genauigkeit von 2—3%) konnen nur aut graphischem Weg 
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gefunden werden. Alle Versuche, direkt auf analytischem Wege eine Lésung zu 
finden, fiihren auf einen unverhaltnismassig grossen rechnerischen Aufwand 
und ergeben dabei oft physikalisch nicht realisierbare Resultate (negative 
Widerstande). Dagegen sind analytische Verfahren geeignet, bereits bekannte, 
gute Naherungsloésungen zu verbessern. 

Alle graphischen Methoden bestehen darin, dass man zwei Kurvenblatter 
aufeinanderlegt und aus der beobachteten Deckungsfigur gewisse Schliisse 
zieht. (Ein durchsichtiger, von unten beleuchteter Ausschnitt auf dem Arbeits- 
tisch erleichtert diese Untersuchungen sehr.) Auf diese Weise kann eine Kurve 
oder sogar eine ganze Schar mit einem Blick erfasst werden, was bei analyti- 
schen Methoden nicht der Fall ist, und ein geiibter Mathematiker wird auf 
Grund seiner Erfahrung schnell die richtigen Schritte einleiten. 

Das einfachste Verfahren ist nun dieses, dass man sich fiir die verschiedenen 
Parameterwerte einen Katalog von Kurven anlegt, aus welchem das am besten 
passende Muster leicht zu ermitteln ist. Da von den Parametern «, B, y aus 
physikalischen Griinden nur ein beschrankter Bereich verwendet werden kann, 
wird der erwahnte Katalog nicht so umfangreich, wie es auf den ersten Blick 
erscheinen mag. Seine Erstellung lohnt sich in solchen Fallen, wo durch die 
gleiche Schaltung viele verschiedene Kurven anzupassen sind. — Andere Ver- 
fahren sind entwickelt worden, die durch gegenseitige Verschiebung von zwei 
aufeinandergelegten Blattern die kontinuierliche Variation von zwei der drei 
Parameter gestatten [8]. Diese Verfahren sind verwandt mit denjenigen, die 
von SvoBopA fiir die Konstruktion von mechanischen Gelenkmechanismen 
angegeben werden [10]. 

Auf diese Art lassen sich leicht Lésungskurven finden, die von den vorge- 
gebenen um nicht mehr als 1—2°% des Maximalwertes abweichen. Die weiteren 
Verbesserungen kénnen nun auf analytischem Wege angebracht werden. Nach 
dem Verfahren, das sich hierzu am besten bewiihrt hat, werden fiir die Para- 
meter nur noch Anderungen zugelassen, die so klein sind, dass ihr Einfluss 
entlang der ganzen Kurve in Form einer Korrektur durch Beriicksichtigung 
der ersten partiellen Ableitungen nach den Parametern erfasst werden kann. — 
Wenn man nun aus der Kurve einige Punkte herausgreift und etwa fordert, 
dass die Summe der Fehlerquadrate an diesen Stellen minimal werde, so erhalt 
man drei lineare Gleichungen fiir die Parameteranderungen, die sich leicht auf- 
ldésen lassen. 

Selbstverstandlich gelten alle diese Bemerkungen nicht nur fiir die Schal- 
tung nach Figur 4, sondern ebenso fiir die vielen anderen méglichen Schal- 
tungen zur Darstellung von Funktionen. 

Es gibt Falle, wo eine gesuchte Darstellung mit linearen Potentiometern 
nicht gefunden werden kann, besonders dann, wenn die vorgegebene Funktion 
im Innern des Bereiches mehr als ein Extremum hat. Haufig fiihrt hier die An- 
bringung von Abgriffen an den Potentiometern zum Ziel. Zwar entfernt man 
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sich hier vom angestrebten Prinzip der Verwendung nur handelsiiblicher Teile ; 
doch lassen sich Abgriffe auch nachtraglich leicht anbringen. Als Beispiel wah- 
len wir die Funktion f(g) = sing, die im Bereich von 0 bis 2 a mit beschrankter 
Genauigkeit darzustellen sei. (Diese Forderung stellt sich zum Beispiel bei der 
Beriicksichtigung des Windes in seinem Einfluss auf ballistische Korrekturen.) 
Figur 6 zeigt eine geeignete Schaltung. Der Kreis kennzeichnet ein Potentio- 
meter, dessen gleichmassig bewickeltes Widerstandselement sich tiber 360° 
erstreckt. Falls die Abgriffpunkte und die Festwiderstainde geeignet gewahlt 
werden, wird die Funktion sing mit einer gréssten Abweichung von 3,5% vom 
Maximalwert dargestellt. Durch Zufiigen von vier weiteren Abgriffen mit ent- 
sprechenden Widerstanden kann der Fehler auf 1% reduziert werden. Dieses 
Verfahren besteht in einer Annaherung der gegebenen Funktion durch einen 
Streckenzug. 


3. Die Erzeugung von Funktionen von zwei Variablen 


Der in Abschnitt 2 formulierte Satz lasst sich auch auf mehrere Variable 
erweitern: Die Schaltung von Figur 3 soll beispielsweise zwei Achsen haben, 
auf denen Potentiometer sitzen und deren Drehwinkel x und y seien. Dann ist 
die resultierende Spannung eine rationale Funktion von x und y. Die Mannig- 
faltigkeit der méglichen Schaltungen ist hier noch bedeutend grdésser als im 
Fall einer einzigen unabhangigen Variablen, und es stellen sich beim Aufsuchen 
einer Anordnung, die eine vorgegebene Funktion realisieren soll, ernstliche 
Schwierigkeiten entgegen, selbst wenn die Existenz einer solchen Schaltung 
erwiesen ist. Eine grosse Vereinfachung entsteht jedoch, wenn man sich auf 
Funktionen beschrankt, die sich als Sechseckgewebe darstellen lassen. Dieser 
Fall sei im folgenden betrachtet. Die dabei zur Verwendung gelangenden Tat- 
sachen aus der Geometrie der Gewebe sind in [2] ausfiihrlich behandelt. 

Die darzustellende Funktion sei z = 2(x, y) und sei durch Niveaulinien in 
einem (x, y)-Koordinatensystem graphisch dargestellt. Die Parallelen zur 
x-Achse seien als x-Schar, die Parallelen zur y-Achse als y-Schar und die 
Niveaulinien als z-Schar bezeichnet. Dieses Bild von drei Kurvenscharen ist 
ein Kurven-3-Gewebe, falls zwei beliebige Kurven von verschiedenen Scharen 
héchstens einen Punkt gemeinsam haben. Wenn sich nun von jeder der drei 
Scharen eine Auswahl von Kurven derart zeichnen lasst, dass an jedem Schnitt- 
punkt genau drei dieser Kurven zusammentreffen, wie dies in Figur 7 gezeigt 
ist, so handelt es sich um ein Sechseckgewebe. Um zu kontrollieren, ob ein 
gegebenes Kurven-3-Gewebe ein Sechseckgewebe sei, ist es jedoch nicht notig, 
im ganzen Gebiet die erwihnte Auswahl von Kurven zu zeichnen; vielmehr 
gentigt das Konstruieren von einigen Brianchon-Sechsecken, von denen in 
Figur 7 zwei verstarkt ausgezogen sind. Von einem beliebigen Anfangspunkt 
aus wird ein Polygonzug gezeichnet, wobei in der aus Figur 7 ersichtlichen 
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Weise entlang einer Kurve der Schar y, z, x, y, 2, x zu fahren ist. Falls man in 
allen Sechsecken, die sich im gegebenen Gebiet zeichnen lassen, wieder an den 
Anfangspunkt zuriickkehrt (wie im oberen Sechseck in Figur 7), so ist das 
Gewebe im betreffenden Gebiet ein Sechseckgewebe. Wenn jedoch der End- 
punkt A’ mit dem Anfangspunkt A nicht zusammenfallt (unteres Sechseck), 
so ist in dieser Umgebung die Sechseckgewebe-Bedingung nicht erfillt. Im 


3 4 5 
6 
/ 
7 
2 
o 
a 
4 
A\ /A’ 
7 Pet 3 4 5 
Fig. 7 
Sechseckgewebe. 


vollstandig gezeichneten Gewebe dussert sich dies so, dass an der betreffenden 
Stelle statt des Schnittpunktes der drei Kurven ein kleines Dreieck entsteht, 
welches bei A deutlich sichtbar ist. 

Wenn nun eine Funktion z(x, y) in einem gewissen Gebiet die Sechseck- 


gewebe-Bedingung erfiillt, so bedeutet das, dass sie sich bei geeigneter Wahl | 


dreier Funktionen X(x), Y(y) und Z(z) wie folgt darstellen lasst: 
Z(2) = X(x) + Y(y). (2) 


Diese wichtige Feststellung bedeutet, dass solche Funktionen durch Verwen- 
dung von Funktionen einer einzigen Variablen dargestellt werden kénnen, wo- 
durch die Schwierigkeiten bei der Verwirklichung mit Potentiometern auf einen 
Bruchteil der urspriinglichen reduziert sind. — Da mit Potentiometern neben der 
Addition auch die Multiplikation leicht ausgefiihrt werden kann, ist es oft 
zweckmassig, Gleichung (2) durch die Substitution 


X(x) = log X"(x), ¥(y) = log ¥'(y), Z(z) = logZ'(y) 


me 
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wie folgt zu schreiben: 
| Z"(z) = X"(x) Y"(y) . 
Hier kann ausserdem ein frei wahlbarer Exponent a eingefiihrt werden: 
[Z'@)}* = [XP LY’ OI. 


Dieser Exponent gestattet es, den darzustellenden Funktionen eine zusatzliche 
Kriimmung nach oben oder nach unten zu verleihen, um sie in eine fiir die 
Anpassung geeignete Form tiberzufihren. 

Alle diese Uberlegungen gelten, wie erwahnt, nur fiir Sechseckgewebe. Es 
zeigt sich nun, dass erstaunlich viele in der Praxis vorkommende Funktionen 
mit zwei Variablen die Sechseckgewebe-Bedingung mit befriedigender Genauig- 
keit erfiillen, obwohl auf analytischem Wege eine Separation nach Gleichung 
(2) in den seltensten Fallen méglich ist. Es lasst sich geradezu aussagen, dass 
von den praktisch vorkommenden Funktionen fast alle diejenigen, die die Be- 
dingung fiir ein 3-Gewebe erfiillen, annahernd auch Sechseckgewebe sind, 
obwohl ein solcher Satz natiirlich nicht allgemeine Giiltigkeit haben kann. Zur 
Auffindung der gesuchten Funktionen X(x), Y(y), 2() muss nach Figur 7 das 
gesamte Sechseckgewebe konstruiert werden. Die Kurven jeder Schar werden 
dann gemass Figur 7 fortlaufend numeriert, und zwar so, dass an jedem 
Schnittpunkt die Nummer der z-Kurve gleich der Summe der Nummern der 
x-Geraden und der y-Geraden ist. Die gesuchten Funktionen sind nun durch 
die Beziehungen zwischen diesen Nummern und den zugehérigen x-, y- und 
z-Werten gegeben. Das Zeichnen des Gewebes bereitet keine Schwierigkeit, 
falls die gegebene Funktion ein exaktes Sechseckgewebe ist. Wenn aber die 
Schliessungsbedingungen nur annahernd erfiillt sind, so muss man bestrebt 
sein, die Ungenauigkeiten (wie bei A, Figur 7) im ganzen Gebiet méglichst 
klein zu halten, was erst nach mehrmaligen Versuchen gelingt und einige 
Ubung erfordert. 

Es ist hier noch eine Bemerkung zum Ausdruck «Ungenauigkeit » notwendig. 
In der Gewebetheorie unterscheidet man Gewebe, die die Bedingung fiir ein 
Sechseckgewebe erfiillen, und solche, die sie nicht erfiillen. « Annihernde» Sechs- 
eckgewebe existieren nicht; denn durch topologische Deformation kénnen Un- 
genauigkeiten wie in Figur 7, wie klein sie auch sein mégen, beliebig vergrossert 
werden, so dass ein solches Gewebe beliebig stark von einem Sechseckgewebe 
abweicht. In unserer Anwendung sind jedoch topologische Transformationen 
nicht ohne weiteres zulassig ; vielmehr bildet die Grdsse der Abweichungen in 
der Gewebestruktur ein genaues Mass fiir die Fehler, die entstehen, wenn die 
betreffende Funktion nach Gleichung (2) dargestellt wird. Wenn diese Fehler 
zum Beispiel 1% des Maximalwertes nicht iiberschreiten, so ist die Feststel- 
lung, das Gewebe sei « annahernd ein Sechseckgewebe», durchaus angangig. — 
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Eine Diskussion einiger Falle von angenaherten Sechseckgeweben sowie eine 
Anregung zu deren Klassifikation ist von SvopopA im Zusammenhang mit 
Gelenkmechanismen gegeben worden [10]. 

Die vorstehend dargelegten Verfahren kénnen nicht angewendet werden, 
wenn die gegebene Funktion z(x, y) von einem Sechseckgewebe stark abweicht, 
ganz besonders dann, wenn nicht einmal die Bedingungen fiir ein 3-Gewebe 
erfiillt sind, indem zwei Kurven verschiedener Scharen mehr als einen Schnitt- 
punkt haben. Dies ist der Fall zum Beispiel bei ballistischen Schusstafeln, die 
neben dem aufsteigenden auch den absteigenden Ast der Flugbahn enthalten. 
Dann ist von Potentiometerschaltungen auszugehen, welche eine entsprechende 
Funktion darstellen, und es muss von vornherein auf eine Separierung nach 
Gleichung (2) verzichtet werden. Um die grossen Schwierigkeiten, die beim. 
Aufsuchen einer geeigneten Anordnung entstehen, etwas zu reduzieren, eignet 
sich folgender Weg: 

Um die gegebene Funktion z(x, y) elektrisch darzustellen, definieren wir 
wieder drei Funktionen X(x), Y(y), Z(z), das heisst, Funktionen mit nur einer 
Variablen, deren elektrische Darstellbarkeit wir als gegeben annehmen und 
daher hier nicht weiter betrachten wollen. Durch Einsetzen in die gegebene 
Funktion z(%, y) entsteht eine neue Funktion Z = Z(X, Y), die als Grund- 
funktion bezeichnet sei. Diese Grundfunktion ist es, die nun dargestellt werden 
muss. Mit einer einzigen Grundfunktion kénnen durch geeignete Wahl von 
X(x), Y(y) und Z(z) sehr viele verschiedene Funktionen z(x, y) dargestellt 
werden. Nun lassen sich die Schaltungen zur Erzeugung von Funktionen mit 
zwei Variablen mit Potentiometern leicht klassifizieren und in ihrem Verhalten 
iibersichtlich darstellen, sofern man sich auf die Verwendung von insgesamt 
zwei Potentiometern beschrankt [11]. Die so erzeugten, verhaltnismassig weni- 
gen Funktionen werden als die zur Verfiigung stehenden Grundfunktionen be- 
trachtet. Wenn nun eine Funktion 2z(x, y) gegeben ist, so sind drei Funktionen 
X(x), Y(y), Z(z) zu suchen, welche z(x, y) in eine der gegebenen Grundfunk- 
tionen tiberfiihren. Der Sinn dieses Verfahrens ist der, eine Beschrankung auf 
wenige, einfache Grundfunktionen zu schaffen und einen méglichst grossen’ 
Teil der Schwierigkeiten auf die Darstellung von X(x), Y(y) und Z(z) iiberzu- 
walzen, welche viel leichter zu beherrschen ist, da es sich um Funktionen mit 
nur einer Variablen handelt. 

Es sei bemerkt, dass die Darstellung von Sechseckgeweben nach Gleichung 
(2) ein Spezialfall dieses Vorgehens ist, indem die einzige Gundfunktion, die 
man zuliasst, die Addition ist. Leider existiert aber bis jetzt kein Verfahren, die 
Funktionen X(x), Y(y) und Z(z) fiir allgemeinere Grundfunktionen in ahnlich 
einfacher Weise zu finden, wie dies im Spezialfall der Sechseckgewebe méglich 
ist. Zum Auffinden einer ersten Naherung ist man daher auf ziemlich umstiand- 
liche, graphische Versuche angewiesen. Verbesserungen kénnen alsdann auf 
analytischem Wege angebracht werden. 
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4. Der Bau von Rechengeraten 


Unter Verwendung der skizzierten Grundschaltungen kénnen Rechengerate 
zar Auswertung bestimmter, einfacher Ausdriicke gebaut werden, welche un- 
bestreitbare , Vorziige gegentiber mechanischen Geraten aufweisen, die mit 
Gelenkwerken, Differentialgetrieben, Kurvenscheiben und Kurvenkorpern ar- 
beiten. Lineare Potentiometer werden in hochprazisen Ausfiihrungen (Genauig- 

_keit bis 1/9) serienmassig hergestellt und sind daher kurzfristig und zu massi- 
gen Kosten erhialtlich. Die gebauten Gerate sind somit schnell und mit be- 
scheidenem Kostenaufwand zusammenstellbar. 

Ist das Resultat durch implizite Lésung einer einzigen Gleichung nur mit 
grossem Aufwand darstellbar, so kann oft durch das Bilden von Zwischen- 
resultaten, das heisst durch das Vorsehen mehrerer Nullabgleiche, eine grosse 
Materialersparnis erzielt werden. Der Vorteil der elektrischen gegentiber der 
mechanischen Darstellung von Funktionen wird besonders gross, wenn in die- 
selbe Gleichung nach Wahl eine von mehreren vorgegebenen Funktionen ein- 
gehen soll, beispielsweise in artilleristischen Rechengeraten fiir die Beriick- 
sichtigung verschiedener Ladungen. Die Umschaltung geschieht in diesem Fall 
rein elektrisch durch Anderung der Festwiderstande, wihrend die Potentio- 

meter dieselben bleiben; im mechanischen Fall ist die Umschaltung von Kur- 
venkérpern mit vielerlei Unannehmlichkeiten verbunden. 

Schliesslich ist noch zu erwahnen, dass solche Rechengerate robust in der 
Handhabung und einfach in der Bedienung sind. Ihr Arbeiten ist unabhangig 
von Spannungsschwankungen der Stromquelle, da alle Resultate durch Null- 
abgleich von briickenartigen Schaltungen gewonnen werden. In vielen Fallen 
geniigt als Quelle sogar eine Trockenbatterie von 4,5 V. 


Der Verfasser méchte Herrn Professor STIEFEL fiir seine Anregung zur Durch- 
fiihrung dieser Untersuchungen sowie fiir seine Unterstiitzung wahrend der 
Arbeiten an dieser Stelle herzlich danken. 
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Summary 


This paper gives some circuits that may be used in computing instruments 
using only resistances and linear potentiometers. Part 1 is a summary of the 
well-known methods of representing the elementary operations. Part 2 is a dis- 
cussion of curve fitting methods for functions of one variable. Part 3 deals with 
curve fitting for functions of two variables; special use is made of the concept of 
grid structure (Sechseckgewebe). Part 4 gives some general remarks regarding 
the computing instruments. 

(Eingegangen: 1. August 1952.) 
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Einige Beobachtungen tiber das Haften von Eis an Oberflaichen 


Von THomMAS BRUNNER!), Weissfluhjoch-Davos?) 


Problemstellung 


Die Untersuchungen sollten einen Aufschluss tiber die Grésse der Adhdsion 
von Eis an verschieden behandelten metallischen und nichtmetallischen Ober- 
flachen ergeben. Ferner sollte abgeklart werden, wie die Oberflachen zu behandeln 
sind, um das Haftvermégen zu vermindern. 


1. Messmethoden 


Die mechanischen Eigenschaften von Eis variieren je nach Herstellungsart — 
ziemlich stark. So wird sich zum Beispiel ein Stiick Eis, das bei — 20° C entstan- 
den ist, anders verhalten, als eines, das langsam bei —3° C gefror. Um vergleich- 
bare Resultate zu bekommen, miissen deshalb die Versuchsbedingungen konstant 
gehalten werden, Da die Resultate aber trotzdem noch stark streuen, fiihren nur 
Reihenversuche zum Ziel. Um verschiedene experimentelle und apparative 
Schwierigkeiten zu umgehen, haben wir fiir die folgenden Untersuchungen den 
statischen Schubversuch gewahlt. 


1) Die Redaktion hat die schmerzliche Pflicht, mitzuteilen, dass der Autor am 1. August 1952, 
am Tage nach der Fertigstellung dieser Arbeit, bei einer Zugsentgleisung tédlich verungliickt ist. 
Herr Dr. M. Dr Quervain, Leiter des Eidgenéssischen Institutes fiir Schnee- und Lawinenfor- 
schung, Weissfluhjoch-Davos, hatte die Freundlichkeit, die Korrekturen der Arbeit zu tibernehmen. 

*) Die Untersuchungen wurden im Auftrage der Schweizerischen Kommission fiir Vereisungs- 
fragen in den Laboratorien des Eidg. Institutes fiir Schnee- und Lawinenforschung, Weissfluhjoch- 
Davos, durchgefiihrt. 
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Gemass den Messungen von FINLAysoN!) an natiirlichem Eis betragt die 
Schubfestigkeit von Eis 7-8 kg*/cm?. Dieser Wert ist niedrig. An sorzfaltic ge- 
reinigten Oberflachen haftet Eis sehr gut. Es wird deshalb eher brechen als sich 
von der Oberflache lésen. Wenn das Haftvermégen von Eis grésser als die Schub- 
festigkeit ist, lasst sich die Adhasion mit diesen mechanischen Methoden prinzi- 
piell nicht bestimmen. 

Hingegen liesse sich mittels radioaktiven Wassers das Haftvermégen einzelner 
H,O-Molekiile messen. Beobachtet man namlich die Desorption der radioaktiven 


Fig. 1 
Aufbau der Scherproben. 


I | 


———-4+_— gesagt 
Yili 
Fig. 2 Fig. 3 
Die Fertigbearbeitung der Proben. Das pordse Die Zugapparatur. Die Scher- 
Eis des mittleren Teiles wird entfernt. spannung wirkt in der Haftebene. 


Molekiile in Funktion der Temperatur und der Zeit (wobei natiirlich fiir die 
Halbwertszeit zu korrigieren ist), so bekommt man leicht die Aktivierungs- 
energie, das heisst das molekulave Haftvermogen von H,O, unabhangig von den 
mechanischen Eigenschaften des Eises. Uber die Aktivierungsenergien von Me- 
tallen auf Metalloberflachen liegen verschiedene mit dieser Methode gewonnene 
Resultate vor); ahnliche Angaben fiir Eis fehlen. 


2. Apparatur 


Auf Metallklétzchen von etwa 2 x 2cm und 1cm Dicke wurde bei zirka 
_4°C destilliertes Wasser aufgefroren. Um diese Klétzchen in die Zugapparatur 
einzuspannen und sie gleichzeitig vor unerwiinschter Warmeiibertragung beim 
Manipulieren zu schiitzen, wurden sie auf Trager aus Cibanit-Isolierplatten auf- 
geschraubt. 

Figur 3 zeigt die Zugapparatur. Die Praparate werden bei B eingespannt. 
Ein Zughaken Z greift unter das Eis E, das auf dem Metallklétzchen M aufge- 
froren ist. Der Zughaken beriihrt die Metalloberflache und wird durch den An- 
schlag A parallel gefiihrt. Der Zug wird mit der Federwaage F ausgetibt. 


1) Zitiert bei: H. N. Dorsey, Properties of Ordinary Water Substance (Am. Chem. Soc. Mo- 


nograph. Ser. 81, New York 1940), S. 428 ff. 
2) H. FRAUENFELDER, Helv. Phys. Acta 28, 347 (1950). 
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3. Priparatherstellung 


Das Wasser wurde im Laboratorium bei —4° C auf die Metalloberflache ge- 
froren, indem ein Klebeband um das Metallklétzchen geklebt und der so entstan- 
dene Raum mit handelsiiblichem destilliertem Wasser von Zimmertemperatur 
aufgefiillt wurde. Um die sehr leicht eintretende Unterkthlung zu vermeiden, 
wurden die Proben nach 60-90 min mit Schnee geimpft und in der Regel 18h 
spater bei gleicher Temperatur untersucht. 

Das Eis wuchs gleichzeitig von der Metalloberflache nach oben und von der 
freien Oberflache nach unten, so dass in der Mitte eine deutliche und immer gleich 
aussehende, durch viele Luftblaschen gekennzeichnete Trennflache entstand. 
Bei Belastung der Probe bricht das Eis meist langs dieser Flache, was uner- 
wiinscht ist. Bei der Mehrzahl der Proben wurde deshalb der obere Teil des 
Eisklétzchens weggesagt. Vgl. Figur 2. Der iibrigbleibende Teil ist glasklar und 
fast frei von Luftblaschen. 


4. Oberflachenbehandlung 


Das Haften von Eis wurde an folgenden Oberflachen untersucht: 

1. Saubeve Oberfldche. Die Praparate wurden mit heissem Wasser, dem ein 
Benetzungsmittel zugesetzt wurde, gereinigt und anschliessend unter fliessendem 
Brunnenwasser geschmirgelt. Dadurch werden die beim Schmirgeln weggehenden 
Partikel fortwahrend weggespiilt. Das die Oberflache benetzende Wasser wurde auf 
die Praparate festgefroren, bevor es verdunsten konnte, wodurch sich ein eventuell 
vorhandener Film nicht auf die Metalloberflache festsetzen konnte. Beim Alu- 
minium bleibt das Oberflachenoxyd sicher auch nach dem Schmirgeln bestehen. 

2. Oxydierte Oberflachen. Die sauberen Oberflachen wurden mit einem 
Schweissbrenner erhitzt. 

3. « Technische » Oberflachen. Die Metalle kamen so zur Verwendung, wie sie 
von der Fabrik geliefert wurden. 

Eine zweite Gruppe von Proben wurde mit einem diinnen Film iiberzogen: 

4, Silikonfett. Auf die saubere Oberflache wurde Silikonfett méglichst diinn 
aufgetragen. 

5. Kunstharz Polyweld, Die sauberen Oberflachen wurden méglichst diinn 
mit Polyweld (IKunstharz auf Polystyrenbasis) bestrichen. 

6. Paraffin. Die sauberen und getrockneten Oberflachen wurden mit handels- 
tiblichem Paraffin tibergossen und nachher mit einem Lappen abgewischt, so dass 
nur ein diinner Film tibrigblieb. 

7. Skilack. Die sauberen Oberflachen wurden mit Skilack bestrichen. 


5. Resultate 


Die Resultate, die an Aluminium, Kupfer und Messing gewonnen wurden, 
werden zusammengefasst, weil sich gezeigt hat, dass die Metallsorte nur von unter- 
geordneter Bedeutung ist. 

a) Bruchverhdltnis 


Da bei allen Proben die Adhasion von der gleichen Gréssenordnung wie die 
Schubfestigkeit ist, bricht immer ein grosser Teil der Proben. Nimmt man an, 
dass die Schubfestigkeit des Eises — es wurde immer unter gleichen Bedingungen 
hergestellt — ungefahr konstant ist, so gibt das Bruchverhaltnis 


ee Anzahl Proben ohne Bruch abgelést 


Gesamtzahl Proben 


einen ersten Anhaltspunkt iiber das Haftvermégen. 
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Aus Tabelle 1 geht hervor, dass bei technischen Oberflachen rund die Halfte, 
bei sauberen Oberflachen jedoch nur zirka ein Sechstel aller Proben ohne Bruch 
von der Oberflache entfernt werden kénnen. An sauberen Oberflachen haftet das 
Eis somit besser als an technischen. Die fiir Ablésung giinstigste Oberflache ist 
eindeutig die mit Silikonfett behandelte. 

Aus Tabelle 1 darf man nicht den falschen Schluss ziehen, dass Skilack, 
Kunstharz und Paraffin ungiinstiger als eine technische Oberflache seien. Die 
nachfolgende Betrachtung zeigt némlich, dass die Schubfestigkeit des Eises keine 
Konstante ist und von der Oberflachenbeschaffenheit abhangt. Wir sprechen 
deshalb lieber von einer «scheinbaren» Schubfestigkeit des Eises. Die Resultate 
in Tabelle 1 lassen sich somit héchstens innerhalb einer Gruppe vergleichen. 


Tabelle 1 
| on one 
: Bruchverhaltnis } Anzahl Proben 
Oberflache B as 
GruppeT 
Sauberh)haru 179%, 46 
Oxydiert(2) . 39% 56 
Technisch (3) . 51% 42 
Gruppe 2 
Skilack (7) . . S594, Sul 
Kunstharz (5) . 399% 33 
Paraffin (6) . . 50% 48 
Silikonfett (4) . 719% 43 


Ve 


b) Schubspannungen 


Nun fragen wir nach den Schubspannungen 7, die ndtig sind, um eine Ober- 
flache eisfrei zu bekommen, unabhangig davon, ob das Eis dabei bricht oder nicht. 
Eine Oberflache bezeichnen wir als eisfrei, wenn mindestens vier Fiinftel der Flache 
vollkommen blank und ohne sichtbare Spuren von Eis sind. Auf dem restlichen 
Fiinftel der Flache diirfen héchstens diinne Scherben, Splitter oder Flecken haften 
bleiben. v bezeichnet den Prozentsatz der Proben, die gebrochen sind, ohne dass 


dabei die Oberflache eisfrei wurde. 


Tabelle 2 
a ok oe eS ee Se Te 
P kg*/cm? Total Proben 
Oberflache - Y x 
Se ee EN en ee rae eee 
Gruppe 7 
Sauber (1) : 9,5 41% 2 
Oxydiert (2) - 8,8 36% a 
Technisch (3) 7,2 15% 


ates LS) Se 


Gruppe 2 
Kunstharz (5) . 6,3 DAs 33 
i 5,9 SIE yh 
Skilack (7) ; 7 3h 
Paraffin (6) .. - 4,6 a> 7 2 
Silikonfett (4) . Shit 2% 


ie Sg ie ea 
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Wieder zeigt sich, dass Eis an sauberen Oberflachen sehr gut haftet. Bei 
9,5 kg*/cm? bleiben immer noch 41% der Proben mit Eis behaftet. Gtinstig ist 
wiederum Silikonfett, wo 3 kg*/cm? Schub geniigen, um praktisch alle Proben 
eisfrei zu bekommen. 


c) Unterschiede zwischen Oberflachen mit und ohne Film 


Zahlt man, wie viele Proben im Schubspannungsintervall 0-1, 1-2, 2-3 usw. 
kg*/cm? ohne zu brechen von der Oberflache abgeschert wurden, so bemerkt 
man, dass sich einerseits die sauberen, oxydierten und technischen Oberflachen 
und anderseits die mit einem Film tiberzogenen Oberflachen unter sich ahnlich 
verhalten. 

Bei den mechanisch behandelten Oberflachen (Gruppe 1) streuen die Resul- 
tate stark, und das Haftvermégen liegt verhaltnismassig hoch. 


Tabelle 3 


Gruppe 1 
mechanisch 
behandelt 
(in kg*/cm?) 


Gruppe 2 
mit Filmen 
(in kg*/cm?) 


Mittleres Haftvermégen. . ..... oa f 4,3 
Mittlere scheinbare Schubfestigkeit . 9,0 4,9 
AnZahLProDeILn 5 a lune ee 155 169 


Bei den mit Film iiberzogenen Oberflachen (Gruppe 2) ist die Streuung klein 
und das Haftvermégen verhaltnismassig gering. Die folgenden Figuren geben eine 


Ubersicht. 


Gruppe 2 96Proben 


20 20 
i- 
15 15 i leru pe 1 98 Proben 
- ss 
a 
c 10 \ 2 70 
A at Gruppe! 57Proben < 
& vhs = 
Cure eh | ‘ Ss 5 
= NU! { N 
r= = 
N 
a 0 0 
= 2 x 
10 15 kg'/cm ) 5 10 15 kg'/ cm’ 
Fig. 4 Fig. 5 
Anzahl Proben, die ohne Bruch abgeschert Anzahl Proben, die gebrochen sind, 
wurden, in Funktion der Schubspannung. in Funktion der Schubspannung. 


Das durch unsere Versuche gemessene Haftvermégen betragt bei mechanisch 
bearbeiteten Oberflachen (Gruppe 1) etwa 8-9 kg*/cm*, bei Oberflachen mit 
Filmen (Gruppe 2) rund die Halfte von diesem Wert, namlich 3,5—4,5 kg*/cm?. 
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Die scheinbare Schubfestigkeit des Eises ist etwa gleich gross und stellt sich somit 
als eine Funktion der Oberflachenbeschafjenheit heraus. — Als wirkliche Schub- 
festigkeit bezeichnen wir diejenige Festigkeit, die man misst, wenn der Bruch 
im Eis und nicht an oder in der Nahe einer Metalloberflache stattfindet. Sie muss 
bei gleichen Versuchsbedingungen eine Konstante sein. 

Die Haftverminderung bei Oberflachen mit Filmen und die scheinbare Schub- 
festigkeit des Eises lassen sich wie folgt erklaren: 

Bei den mechanisch behandelten Oberflachen der Gruppe 1 gibt es eine 
grosse Anzahl von isolierten Zacken, Spitzen, «Graben» und «Talern», das heisst, 
die wirksame Oberflache ist bedeutend grésser als die scheinbare. Wird nun von 
aussen eine Schubspannung angelegt, so verteilt sie sich auf eine grosse Flache. 
Die daussere Schubspannung kann also ziemlich gross werden, bis die kritische 
mikroskopische Schubspannung iiberschritten wird. Sobald dies irgendwo ein- 
tritt, lést sich dort die Bindung, oder es bricht dort das Eis. Es haftet deshalb 
an einer verkleinerten Fliche, und die Spannung nimmt zu, ohne dass sie von 
aussen vermehrt worden wire. Dadurch lost sich oder bricht das Eis an weiteren 
Stellen usw. und fallt schliesslich von der Oberflache ab. Den gleichen Vorgang 
hat BowpEN?) beim Ubergang von der Haftreibung zur Gleitreibung bei Metallen 
festgestellt. 

Wenn die Oberflachen aber mit einem Film iiberzogen werden, wird ein Teil 
der «Graben» und «Taler» im Metall aufgefiillt, das heisst die effektive Oberfldache 
vervingert. Bei zunehmender dusserer Spannung wird bei der verminderten effek- 
tiven Oberfliche die zuldssige mikroskopische Schubfestigkeit und das Haft- 
vermégen rascher erreicht; scheinbare Schubfestigkeit und Haftvermogen er- 
scheinen deshalb geringer. 

Es leuchtet ein, dass die effektive Oberflache mittels Filmen héchstens auf die 
makroskopische, messbare Oberflache verringert werden kann. Methoden, die mur 
die effektive Oberflache verringern, sind also nur von begrenzter Wirksamkeit. 


d) Das Verhalten der einzelnen Metalle 


Das Haftvermégen von Eis kann durch sorgfaltige Reinigung der Oberflache 
vergréssert werden. Beim Kupfer hat das Schmirgeln den gréssten, bei Aluminium 
den geringsten Einfluss. Das riihrt wahrscheinlich daher, dass durch das Schmir- 
geln beim Kupfer und zum Teil auch beim Messing die Oberflachenoxydschicht 
entfernt wird, wahrend sie beim Aluminium bestehen bleibt. 

Im iibrigen verhalten sich die Metalle nicht einheitlich. Zwanzig Proben pro 
Oberflachenart und Metallsorte reichen wegen der grossen statistischen Schwan- 
kungen nicht aus, um eine sichere Aussage iiber das Haftvermégen zu geben. 
Immerhin kann mit Bestimmtheit gesagt werden, dass von den Metallen Kupfer, 
Aluminium, Messing keines speziell giinstig oder ungiinstig ist. 


6. Andere Methoden zur Verminderung des Haftvermogens 


In der Literatur sind die Hinweise aut die Verminderung des Haftvermogens 
zwar zahlreich, aber meistens wenig aufschlussreich. Fast alle Autoren beschran- 
ken sich darauf, die Methoden (Fette, Ole, Alkohole, Pasten und Salze) und ihre 
Vor- und Nachteile aufzuzahlen. / 

Alkohole, tiberhaupt Fliissigkeiten, ergeben zwar eine geringe Haftung, hin- 
gegen werden zur E/isfreihaltung grosse (und zum Beispiel bei Flugzeugen schwer 
zu transportierende) Mengen benotigt. 


or ., F. P. BowvEN, Friction, Nature 166, 330 (1950). 


ZAMP III/30 
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Fette und Pasten sind giinstiger, da der mengenmassige Verbrauch geringer ist. 
Hingegen werden sie durch den Regen rasch ausgelaugt und miissen in kurzen 
Intervallen erneuert werden, was verhadltnismadssig umstandlich ist. 

Salze. SmMitH!) hat festgestellt, dass die Adhdsion von Eis auf mit Wachs be- 
handelten Aluminiumoberflachen betrachtlich verringert werden kann, wenn das 
Eis aus einer schwachen Lésung gewisser Salze gefroren wird. So vermindert zum 
Beispiel eine 1/1000-normale NaCl-Lésung die Haftung um 50%, eine solche von 
NaNO, um 85% und von Th(NO,), um 97%. Gelingt es, solche Salze in die zu 
schtitzenden Oberflachen so einzulagern, dass sie nicht zu rasch ausgelaugt werden, 
so wird die Adhasion wirksam vermindert. SmitH hat dieses Problem mindestens 
teilweise dadurch gelést, dass er die Salze mit Formvar gemischt hat und die Form- 
varfolien auf die Oberflache klebte. Diese Methode scheint momentan die beste 
zu sein; es ist uns jedoch nicht bekannt, wie sie seit 1946 weiterentwickelt wurde. 


Schlussbemerkungen 


Verschiedene interessante Fragen wurden von uns nicht abgeklart. So zum 
Beispiel die Abhangigkeit des Haftvermégens von der Temperatur. Nach LouGu- 
BOROUGH?) soll letzteres eine lineare Funktion der Temperatur sein. — Auch die 
Art der Eisbildung wird von Einfluss sein. Raschgefrorenes Eis ist feinkristallin 
und haftet wahrscheinlich schwacher als langsam gefrorenes. Ebenfalls sollte 
die Abhangigkeit der Resultate von der Priifmethode untersucht werden. Sicher 
geben dynamische Methoden andere Resultate als statische. 


Summary 


The adhesion of ice grown on different surfaces is investigated by means of 
a static shearing test. 300 samples all made at similar conditions are tested. The 
adhesion on clean metal surfaces is found to be 8 kg*/cm?; it is therefore of about 
the same magnitude as the shearing stress of ice. A surface-film of silicone grease 
reduces the adhesion to one half. The apparent shearing stress of ice is also 
diminished by the same amount. Both effects are attributed to a reduction of 
the true surface. 


(Eingegangen: 4. August 1952.) 


Weitere Bemerkung zu iy eb(x+acosx) dx 


Von Peter Henricr, Washington, D. C.3) 


In einer kiirzlich in dieser Zeitschrift erschienenen Note‘) gibt K. EmpEn 
eine Lésung fiir das im Titel genannte unbestimmte Integral in Gestalt des 
Produktes einer Exponentialfunktion und einer Fourier-Reihe, deren Koeffizien- 
ten er rekursiv durch Auflésung eines unendlichen Gleichungssystems berechnet. 


1) R. SmiruH-JOHANNSEN, AAF Technical Report 5539, General Electric Co. U.S.A. (1946). 

2) D. L. Loucusorovucn, Reduction of the Adhesion of Ice, J. aeron. Sci. 13, 126 (1946). 

8) American University. Die folgenden Untersuchungen wurden gemass einem Vertrage zwischen 
der American University und dem National Bureau of Standards, Washington, ausgefiihrt. 

*) K. Empen, Eine Lésung fiir [ e®(*+4 008) dx, ZAMP 2, 289-292 (1951 ) 
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H. RUTISHAUSER!) schreibt in einer erganzenden Notiz diese Fourier-Reihe 
komplex und gibt fiir die Koeffizienten eine Kettenbruchentwicklung an. Es ist 
indessen durch Wahl einer andern Variabeln méglich, das fragliche Integral in 
einfacher Weise in eine nach Produkten Besselscher und trigonometrischer 
Funktionen fortschreitende Reihe zu entwickeln, deren Konvergenzeigenschaften 
sich im Gegensatz zu denjenigen der von den genannten Autoren angegebenen 
Reihen leicht tiberblicken lassen. 

Wir bendtigen zur Herleitung dieser Reihe das folgende bekannte Resultat 
aus der Theorie der Besselschen Funktionen: Ist a €(—1, 1) und ¥ eine reelle 
Variable, so wird die zu 


= 4—asinx (1) 


inverse Funktion dargestellt durch die samt allen ihren Ableitungen absolut und 
gleichmassig konvergente Reihe 


zat SH hlna) sin (” y) , (2) 


n=1 


wo J,(ma) die n-te Besselsche Funktion bedeutet?). Damit ergibt sich fur 
a€(—1, 1), fiir komplexes b mit Rb > 0 (was keine wesentliche Einschrankung der 
-Allgemeinheit bedeutet) und fur beliebiges reelles X, wenn wir das zwischen den 
Grenzen —co und X genommene Integral mit E(a, b; X) bezeichnen, 


a8 
E(a, b; X) = [ eo(v+acosx) dy 
eo) 
X—n/2 
(zt [ eb [x+m/2+a cos (x +/2)] ax 
5 
X 2/2 
= eub/2 i e0(¥—asin x) ax 
We) 
X —n/2+acosX 
ax 
— pmb/2 by d 
= en / /[ é dy y 
—oo a 
X —n/2+acosX foe) 
= enbj2 fh eby a 2 pa J,(m @) cosn y} dy 
— oo n=1 


und daher, wenn gliedweise integriert und Y = X — n/2 + acosX gesetzt wird, 


be nsinn Y + bcosn Y 
Birk eae ean [p42 D, Jnl a) ets |. (3) 


1) H. RUTISHAUSER, Bemerkungen zur Arbett von K. Emden etc., ae *, re oe 

2) G. N. WaTSON, A Treatise on the Theory of Bessel Functions, 2. Aufl. Seago 4 pe cee 
bridge 1944), S. 552. Gerade durch das Problem, fur die zu (1) inverse — einen ces Bhs en. 
Ausdruck zu finden, wurde BEsseL auf die spater nach ihm benannten Funxtionen gern’. 
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Da fiir grosse 2 ‘i 
e-n (artgh Vi-a— V1—a?) 3 


Vanni ee 


ist die Konvergenz von (3) besser als diejenige einer geometrischen Reihe mit 
dem Argument exp (— artgh Vi —at+ V1 — a?). 


Jn(n a) ~ 


Summary 


Be 
For the function E(a,b;X)= f[ e(*+@c0s) dy (Rb> 0) a series representation 
—0o 
in term of products of Bessel and trigonometric functions is given, valid for real X 
and real a € (—1, 1). 


(Eingegangen: 8. August 1952.) 


On the Existence of Certain Solutions of a 
Nonlinear Differential Equation 


By GrEorGE SEIFERT, Lincoln, Nebraska, U.S.A. ?) 


The equation 


do : 
eee Skee a 1 
ie + a Psind = 8, (1) 
where « and f are positive constants, is of interest in the study of the oscillations 
of a synchronous motor. If we put 

do dy : 

eee then ae B—sinOd—ay, 


and dividing the latter equation by the former, we obtain the equation of the 
phase trajectories 
dy  B—sinOd—ay A 
do yo i (2) 
It is known) that for fB = 1 there exists a solution y(0) = 0 of equation (2) 
such that (8) = (0 + 22) for all 6. The solution 6(t) of equation (1) correspond- 
ing to such a 7(0) is called a periodic solution of the second kind. If 0< fii ake. 
there exists an a) = (8) such that the conditions «Xa and a> % imply 
respectively the existence and nonexistence of periodic solutions of the second 
kind of equation (1). F. TRIcomI’), who has shown the existence of such an a, 
has given bounds on a, and L, AMERIO4), in a study of the stability of solutions 


1) G.N. Watson, l.c., S. 243, 
*) University of Nebraska. 
3) F. Tricomt, Integraszione di un’e 
norm. sup. Pisa [2] 2, 1-20 (1933). 
4) L. Amerio, Determinazione delle condizione di stabilita per 
sante l’elettrotecnica, Ann. Mat. pura appl. [4] 30, 75-90 (1949). 


quaztone differenziale presentatasi in elettrotecnica, Ann. R. Sci, 


gli integrali di un’ equazione interes- 


a 
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of equation (1), has given an upper bound on «) which is lower than that of 
TRICOMI’S. 
It is the purpose of this note to derive another upper bound on a, and to 


show that it is smaller than the one given by AmeErio. To this end, we first 
prove the 


Theorem. 1f0S B < 1, and if 


(3 7 — 4 69) B + 2 (1 — cos 4) la 


x> au = | Zh) 


. 


where sin Ch , 0 < 6,< 2/2, then equation (2) can have no nonnegative solu- 
tions ¥(0) of period 2 z in 9, and all solutions @(¢) of equation (1) tend to a finite 
limit as t > + oo. 


Proof. Assume that there exists a solution 7(@) = 0 of equation (2) such that 
y(0) = y(0+ 22) for all 6. Then by equation (2), ¥(6) must satisfy the equation 


7 : sn [v2(0)] = B — sin@ — a y(6). (3) 


Substituting y = y(@), integrating this from any fixed 0 to 0 + 22, and dividing 
by «, we obtain 


a 


O+27 

Y) 
if; CE eee 
6 


From this, it follows that the minimum value m of y(@) must be such that #< B/a. 
We also note that (0) = m implies 


6=6,+2an, n=0, +1, £2,.- where zn — 06,5 0,,< 2. 


This is a consequence of the observations that the locus of points in the (9, ¥) 
cartesian plane for which dy/d6 = 0 is the curve 


> pS sine 
JS s ’ 


and that the region included between this curve and the 6 axis is the region in 
which all phase trajectories have positive slope. Hence if B = 0 we have a contra- 
diction and conclude the theorem. We note, incidentally, that in this case ay = 0, 
which is the best possible value. Accordingly, in what follows we may assume 


pe 0. ; 

Let us consider the solution y(0) of equation (2) for which y(z) = B/«. Since 
y(O) > B/a, we have, clearly, that (4m) > (9m): and consequently, y(@) > 7(9) 
for all @. Define the function /(@) as follows: 


peste Pas » 

od 

Ue oy (Oe Oy) t0m 2m ee 0= 855 
a (0, — 0) for : SOS 4% 

where 0, = 7— %, A= OF 2x2, and 6,= 0,+22. Since 


p.m 46) BE2(L— 0080) < gs 
0, — 9 qt (m6 = 2 Oo) 


trethave thatsx, (0) — 0) => Bia hence, /(0,) 28/a, 


? 
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Assume for a moment that (6) < {(@.). We will show that in this case there 
exists a &, 0. S& < 63, for which y(&) = 0. For if not, the conditions y(@) > 0 
and y’(0)< 0 for 6. 0X 6, would imply the intersection of the trajectory 
y = y(0) with the line y = « (8; — 6) at a point 0= n, 0.< 7 < 43, such that 
—«% S y(n). But 


SA -aaen 10 FO ORM, <a oe (4) 
and we have a contradiction. Hence, the assumption (4) < f(@.) implies the 
existence of a § for which y(£) = 0. But since y(@) > ¥(@) = O for all 0, we con- 
clude that 4(6,) < /(0,) is an impossibility; hence we must have y(O,) = f(O.). 

We show next that y(0)=f(0) for 1<6@< 5. Since y(x) = B/a and~ 
v(62) = f(9:), we have, clearly, that »(6) > /(@) for ~<@06,. Now if for 
0,365 6, we do not have y(6)=>«(6;— 6), then either there exists a é, 
6.< €< 6;, such that y(&) = 0, or there exists an n, 0,< < 0, for which 
y’(n) = —«. Again (4) rules out the latter possibility, while the argument of the 
previous paragraph rules out the existence of the £ We have, then that 
y(9) = f(6) for ~< 0< 4. 

Now by equation (2), y(@) must satisfy equation (3), and integrating from a 
to 0, we obtain 


ee — x) + cos cua (6) do + 2-7) 
> = a9 ; y aa ia 


< £(0@,+ 2) +1-—cos@ 


~ a [Za + 0 (0 ~ 0) (0-22) + * (0,— 09°) + 2 (,— 0 


w 


= $ (32-46) + 1—costy— * a (a— 26). 
Hence, by the hypothesis of the theorem, we obtain y*(03)< 0, and we see that 
the assumption of a y(0) => 0 such that y(0) = 9(8 + 22) for all @ leads to a 
contradiction. 
That the nonexistence of such periodic solutions 7(@) of equation (2) implies 
that for any solution 0(t) of equation (1) either 
lim 6(t) = 0,+22an or lim 6() = 6, + 22m 
t—0o t—>0o 
for some integer m or ” is a consequence of the work of L. AMERIO4), This com- 
pletes the proof of the theorem. 
We show, finally, that AMERIO’s upper bound on a»; namely 


(B + 1) [4 cos6, — 2 (% — 2 09) B}-1/2, 


is greater than ay. If we put y = a/2 — 05, we have, clearly, only to show that 
for 0< p< 2/2, 


I 
2 (cos m + 1*ze> [(% + 4) cosp+ 2(1— sin ¢)} [sing — pcosq], 


1) L. Aerio, Determinaszione delle condizione di stabilita per gli integrali di un’ equazione inter- 
essante Velettrotecnica, Ann. Mat. pura appl. [4] 30, 75-90 (1949), 
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which we establish proving that for 0< gy < a/2, 
rt 
OP 2i(¢) = 4 7 singp + pcosp>0, and 
2 


(ii) g2(v) = 2 (cosp + 1)? — (a9 + 4) cosm— 2 (1— sind)>0. 


To prove (i) we simply note that g,(0)=0, g,(a/2)> 0, g{(0)>0 and 
£1(~) < 0 for 0< » < a/2. To prove (ii), we first note that since sin? y < sin y, we 
have that g,(y) = (4— 1— 49) cosp+ 2. Consequently, g,(0)> 0, g,(a/2)> 0, 
and we note that at » = m, where 
4—2 

at 


COt Py = Po — 


the function g,(g) assumes its minimum value in 0< g< a/2. It is easily estab- 
lished that gy is approximately 0:94, and that since g.(qyp) > 0, the proof of (ii) 
is complete. 


Zusammenfassung 
Die Differentialgleichung 
ad*d d : 
ap Lo se Psing = Bp, 


wo « und £ positive Gréssen sind, kommt in der Elektrotechnik bei der Beschrei- 
bung der Schwingungen einer Synchronmaschine vor. Setzt man y = dd/dt, so 
erhalt man dy/dt = Bf — sind — « y, und aus diesen zwei Gleichungen folgt 


db y 
Man beweist: Ist B< 1 und a> ay , wo 


(3a — 4 8) B+ 2 (1 — cosB)) 
<a me (% — 2 Ho) ‘ 
sind, so existieren keine Lésungen ¥(#) der vorigen Differentialgleichung mit 


den Eigenschaften 
(I) y(8) 0, und (Il) 74) = (0 + 22) 


dy B-—sin®—ay 


4 


sind, = 8B, 0<%< 5 


fiir jedes #. Ferner wird gezeigt, dass der Wert a kleiner ist als der von F. TRI- 
comi und L. AMERIO angegebene. 
(Received: July 14, 1952.) 
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8. Internationaler Kongress fiir theoretische und angewandte Mechanik 

Der Kongress fand vom 20. bis 28. August 1952 an der Universitat Istanbul 

statt. Er ist von einem lokalen Organisationskomitee unter der Leitung von 


Prof. Kertm Eri in vorbildlicher Weise vorbereitet und durchgefiihrt worden 
und vereinigte gegen 600 Teilnehmer aus fast allen freien Landern der Erde. 


472 Varia — Miscellaneous — Divers ZAMP * 


Folgende Hauptvortrage wurden gehalten: R. Courant, Reciprocal Varia- 
tional Problems and Applications to Problems of Equilibrium; J. C. HUNSAKER, 
Social Aspects of Aeronautics; A. KANTROWITZ, Physical Phenomena associated 
with strong Shock Waves; A. LICHNEROWICZ, Méthodes tensorielles en mécanique; 
J. von NEumMANN, Automatic Computation and Hydrodynamics; Sir G. TAYLOR, 
Hydrodynamic Theory of Detonating Explosives; W. TOLLMIEN, Uber Schwingungen 
in laminaren Strémungen und die Theorie der Turbulenz. 

Rund 350 Referate verteilten sich auf folgende fiinf Sektionen (von denen die 
beiden ersten weiter unterteilt werden mussten) : 

I. Elastizitat, Plastizitat und Rheologie (Sektionsvortrage: H. GEIRINGER, 
Some Recent Developments in the Theory of an Ideal Plastic Body; W. PRAGER, 
The General Theory of Limit Design; M. REINER, Rheology of Second Order © 
Effects in Elasticity and Hydrodynamics) ; 

II. Hydromechanik (M. BurGERs, Some Statistical Problems connected with — 
Solutions of a simple Non-lineay Equation; F. N. FRENKIEL, Turbulent Diffusion; 
G. TEMPLE, The Dominant Derivative Method in Aerodynamics) ; 

III. Mechanik fester Kérper (A. SIGNORINI, Quelques recherches récentes de 
dynamique des solides) ; 

IV. Statistische Mechanik, Thermodynamik, Warmeleitung (M. Roy, Du 
combustible a la propulsion) ; 

V. Mathematische Hilfsmittel und Rechenverfahren (N. Mrnorsky, Sur la 
méthode stroboscopique et ses applications; M. Piconr, Sur une méthode pour le 
calcul des valeurs propres et des fonctions propres et sur son application & l’equation 
de Schrodinger). 

Der nachste Kongress ist fiir 1956 in Briissel vorgesehen. Hans Ziegler 


Herbsttagung der Schweizerischen Physikalischen Gesellschaft 
vom 24. August 1952 in Bern 


Berichte tiber angewandte Physik und Mathematik 


Uber die Temperatur verteilung an angestrémten Zylindern, von J. AcKr- 
REG, Par Zurich. 

Wie man aus Versuchen von ECKERT, WEISE und von RYAN weiss, ist die 
Temperaturverteilung an einem parallel angestrémten, warmeisolierenden Kreis- 
zylinder sehr unsymmetrisch, so dass auf der Riickseite im Ablésegebiet relativ 
ttefe Temperaturen sich vorfinden. Es wird der Versuch gemacht, diese unerwar- 
tete Absenkung auf der Grundlage einer nichtstationaren Potentialstrémung mit 
Einzelwirbeln im Abléseraum zu verstehen. Fiir die Temperatur-Mittelwerte, 
die bisher allein gemessen wurden, ergeben sich sehr einfache Beziehungen, die 
bei verniinftiger Wahl der Zirkulation und mit Einfiihrung des nicht vollstandi- 
gen Temperaturriickgewinnes in Grenzschichten mit den Messungen recht befrie- 
digend iibereinstimmen. Einige Besonderheiten der Temperaturen im Nachlauf 
und in seiner unmittelbaren Umgebung sind in Ubereinstimmung mit den 
Ryanschen Messungen. 


Ein ausfiibrlicher Bericht wird voraussichtlich in der ZAMP erscheinen. 
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The Secondary Flow in Curved Pipes, von R. W. DertrRa, Ziirich}). 


Fluid, with a nonuniform velocity distribution, flowing in curved pipes 
produces a secondary fluid motion in the cross-section of the pipe perpendicular 
to the principal flow. Under the action of the centrifugal force, the fluid particles 
with the higher axial velocity are displaced toward the outside of the bend. To 
satisfy the continuity conditions within the cross-section of the pipe, the slower 
moving fluid particles, which consequently are acted upon by a smaller centri- 
fugal force, are displaced toward the inside of the bend. This resulting secondary 
flow distorts the axial velocity distribution and in itself represents a loss of 
energy from the principal pipe flow. 

A study of this phenomenon in curved constant-area pipes of large radius of 
curvature and small angle of bend is made. Using the inviscid-fluid approxima- 
tions and considering the secondary flow as a small perturbation superimposed 
on a principal flow, a theory for predicting this secondary flow is developed. By 
considering the centrifugal force as a body force acting on each fluid particle, its 
‘magnitude being proportional to the particles velocity squared, a single ordinary 
differential equation, satisfied by the secondary flow stream function, is obtained 
from the equations of motion. In this differential equation the secondary flow 
exhibits a strong dependence on the inlet velocity distribution. A solution for the 
secondary flow stream function is found as a finite series. Superposition of this 
secondary motion on the principal flow produces a distortion of the axial velocity 
distribution as a function of the angle of bend. 

A comparison of the results of this theory with experiments conducted at the 
Institut fiir Aerodynamik, E.T.H., Zurich, shows good agreement. The loss due 
to this secondary flow, while small for small angles, increases as the square of 
the angle of bend. Loss calculations for an elliptical cross-sectioned pipe, whose 
ratio of major axis to minor axis is 2:1, indicate the loss for the pipe, oriented 
with its major axis parallel to the radius of curvature of the bend, is of the order 
of 5 times as great as the pipe oriented with its minor axis parallel to the radius 
of curvature. 


Schubvermehrung durch Strahlmischung, von Z. PLaskowskI, Ziirich’). 
Die aerodynamische Schubvermehrung durch Strahlmischung wird zunachst 
theoretisch erfasst. Die Berechnungsresultate, die fiir inkompressible Stromun- 
gen unter Vernachlassigung der Wandreibung und fiir vollstandige Mischung im 
Endquerschnitt des Schubvermehrers ermittelt wurden, zeigen eine starke Ab- 
nahme der Zusatzschiibe mit zunehmender Fluggeschwindigkeit. Interessante 
Kraftverhdltnisse liegen theoretisch am Stand fiir relativ zur Treibdtise grosse 
Schubvermehrer vor. y 

Inkompressible Standversuche, in den Jahren 1948/49 im Institut fiir Aero- 
dynamik ETH. durchgefiihrt, zeigten, dass fiir einen Schubgewinn in der Gros- 
senordnung von 25 bis 40% grosse Mischrohrlangen benétigt wurden, die einen 
praktischen Einsatz des Schubvermehrers erschweren. 

Die Versuche ergaben starke Abhangigkeit der Schubvermehrung von der 
Vollstandigkeit der StrahIlmischung, die durch Mischrohrlangen und Querschnitts- 
verhdltnisse beeinflusst wird. Man steilte ebenfalls einen Einfluss der Form des 
Mischrohreinlaufes fest, sobald es sich um Vermeidung der Ablosung und Aus- 
bildung der fiir die Schubvermehrung massgebenden Saugkraft handelte. 


1) Institut fiir Aerodynamik, ETH. 
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Durch Unterteilung des Treibstrahls in vier kleine Diisen und bei passender 


Diisenverteilung konnte ein Fortschritt gegeniiber einem einfachen Treibstrahl 


erreicht werden. Es wurde eine Abhangigkeit der optimalen Diisenverteilung vom 
Querschnittsverhaltnis von Treibdiisen und Mischrohr festgestellt. 

Der héchstgemessene Wert des Schubgewinnes betrug 40 % bei einer extremen 
Anordnung. 

Ein ausfiihrlicher Bericht wird vermutlich in der ZAMP erscheinen. 


Uber thermische Effekte bei Resonanzrohren, von H. SPRENGER, Ziirich}). 


Wenn der Gasinhalt eines relativ starren Hohlraumes periodisch komprimiert 
und expandiert wird, kénnen irreversible Vorgange Temperatursteigerungen ver- 
ursachen, die ein Vielfaches der adiabatischen Temperaturerh6hung betragen. 
Dieser Effekt tritt besonders ausgepragt in Erscheinung, wenn das Gasvolumen 
Resonanzschwingungen mit grossen Amplituden ausfiihrt. Fiir die Untersuchun- 
gen sind Rohrresonatoren verwendet worden, die sich mit Gasstrahlen leicht an- 
regen lassen. Es zeigte sich, dass mit schwachen instationéren St6rungen im 
Gasstrahl, intensive Schwingungen sich auch bei unterkritischen Diisendruckver- 
haltnissen erzeugen lassen — beispielsweise durch die Turbulenz eines quer iiber 
die Diisenaustrittsflache gespannten feinen Fadens. Mit dieser Massnahme kann 
sogar die Wirkung von Gasstrahlen, die bei iiberkritischen Druckverhaltnissen 
expandieren, verstarkt werden — eine Verbesserung, die sich vielleicht auch bei 
den Ultraschallgeneratoren von J. HARTMANN anwenden liasst. Druck-, Tempe- 


ratur- und Schallstarkemessungen sind an zahlreichen Resonanzrohren im Fre- 


quenzbereich 50 bis 3000 Hz mit und in Luft ausgefiihrt worden unter Verwen- 
dung konvergenter Diisenformen mit Druckverhiltnissen bis 7:1 und mit einer 
Laval-Diise fiir die Anstr6mmachzahl 1,7. Qualitative Versuche zeigten, dass die 
Erwaérmungen abgeschwacht auch bei Ultraschallfrequenzen auftreten kénnen. 
Eine optimale Temperatursteigerung von 450° C wurde gemessen am verschlos- 
senen Ende eines diinnwandigen Metallréhrchens von 0,3 cm Durchmesser und 
10 cm Lange (900 Hz) bei koaxialer Anstr6mung aus 1,2 cm Entfernung durch 
eine Diise von 0,3 cm Durchmesser. Es wurde dabei Pressluft von Raumtempera- 
tur und 4 b Uberdruck beniitzt. Durch Isolieren des Resonators und Vorwarmen 
des Druckgases lasst sich die Temperaturerhdhung weiter steigern. Mit den 
thermischen Effekten in Resonanzrohren kénnen exotherme Prozesse eingeleitet 
werden. Zur Demonstration des Vorganges geniigt schon ein Stiick Holz mit 
einer Anbohrung. Wenn gegen diese Offnung ein Luftstrahl (3 bis 4 atii Druck) so 
gerichtet wird, dass intensive Resonanzténe zu héren sind, verbrennt das Holz 
von innen heraus unter lebhafter Rauchentwicklung und Funkenbildung. 

_ Wird in der Nahe des geschlossenen Endes eines Resonanzrohres eine kleine 
Offnung angebracht, so dass ein Teil des sich im Hohlraum ansammelnden Warm- 
gases entweichen kann, so erfahrt der aus dem Rohr riickfliessende Treibgasstrahl 
eine der abgegebenen Kompressionsarbeit entsprechende Abkiihlung unter die 
Drosseltemperatur. Resonanzrohre kénnen wie Ranquesche Wirbelrohre benutzt 
werden zum Trennen expandierender Gase in erwarmte und abgekiihlte Teil- 
stréme. Die Charakteristiken eines «Wirbelrohres ohne Blende» und eines Reso- 
nanzrohres einfachster Bauart sind unter gleichen Betriebsbedingungen gemessen 
worden. Die Ahnlichkeit der Kurven deutet hin auf eine physikalische Ver- 
wandtschaft der beiden Separationsvorginge. 


1) Institut fiir Aerodynamik, ETH. 


— 
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Uber das radiale Gleichgewicht in axialen Turbomaschinen, von P. pE 
HALLER, Winterthur?). 


Wenn in einer axialen Turbomaschine die Verteilung der Umfangskomponente 
der absoluten Geschwindigkeit nicht dem Drallgesetz c,7v = const folgt, kann 
die Stro6mung nicht mehr parallel zur Achse verlaufen, da Druck- und Massen- 
krafte nicht im Gleichgewicht stehen. Es wird an Hand des Beispieles der Stré- 
mung hinter einem Leitrad mit konstantem Austrittswinkel gezeigt, dass in den 
fiir die Praxis wichtigen Fallen die radiale Komponente der Geschwindigkeit 
klein bleibt und dass ein neues Gleichgewicht mit veranderten Verteilungen der 
Axial- und Umfangskomponenten sehr rasch, nach ungefaéhr einer Kanalhéhe, 
erreicht wird. Gleichzeitig werden einige Eigenschaften der Differentialgleichung 
und ihrer Lésungen diskutiert. 


Knickung gerader Stabe durch Druck und konservative Torsion, von 
Max BEck, Ziirich?). 


Deformiert sich der durch das Torsionsmoment W belastete Stab (Figur 1), 
so ist zu erwarten, dass der Momentvektor seine urspriinglich axiale Richtung 
andert. Beim Ausknicken des Stabes sind die Deformationen zuerst klein, und 
die Beziehungen zwischen W und dem Tangenten-Einheitsvektor i(¥’, y’, 1) kén- 
nen allgemein durch den linearen Ansatz 


Woe W (Ax = By), Woe WC # — Dy’); Wy= W 


wiedergegeben werden. Fiir die mathematische Behandlung eines Knickproblems 
ist es wichtig, zu wissen, ob die Belastung konservativ sei oder nicht’). Die Be- 
rechnung der virtuellen Arbeit ergibt, dass das Torsionsmoment konservativ ist, 
wenn die Bedingung A+ D =1 erfiillt ist. Die einfachsten konservativen Mo- 
mente erhalt man fiir d = 1; B=C=D=0: 


= 


W = W (#’, 0,1) : quasitangentiales Moment , 


und fiir Ad = D=1/2; B=C=0 


W=W (5 ae : ae 1) : semitangentiales Moment. 


Dagegen sind das axiale Moment (A = B=C = D = 0) und das tangentiale Mo- 
ment (A = D=1; B = C = 0) nicht konservativ. . 

Um zu erfahren, wie sich der Momentvektor bei der Deformation verhalt, muss 
man auf die das Moment erzeugenden Einzelkrafte zuriickgehen. H. ZIEGLER *) 
hat gezeigt, dass man mit der Annahme, dass die Krafte ihre urspriingliche Rich- 
tung beibehalten, konservative Momente erhdalt. Die einfachsten und wohl wich- 
tigsten Anordnungen der Krafte fiihren auf quasi-, semi- oder pseudotangentiale 
Momente. H. Z1EGLER findet in allen drei Fallen fiir alle Arten der Lagerung 
plausible Knickmomente, wahrend die bisher tbliche, wohl selten verwirklichte 
Annahme, der Momentvektor bleibe beim Ausknicken axial, fiir gewisse Lage- 
rungen auf unbefriedigende Resultate fiihrt. 


1) Gebriider Sulzer AG. ; or 

2) Assistenz fiir technische Mechanik, ETH. Riek 
5 H. Z1EctErR, Die Stabilitdtskriterien der Elastomechantk, Ing.-Arch. 20, 49 (1952). 
4) H. ZimcierR, Knickung gerader Stibe unter Torsion, ZAMP 3, 96 (1952). 


' 
- 
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Diese Rechnungen sind nun auf den auf Druck und Torsion beanspruchten 
Stab sowie auf den Stab mit zwei verschiedenen Biegesteifigkeiten und den ver- 
wundenen Stab ausgedehnt worden. Von den Resultaten sind hier nur die Knick- 
lasten fiir den durch eine axiale Druckkraft und ein semitangentiales Moment 
beanspruchten Stab wiedergegeben (Figur 2). Ein Stab ist stabil, solange seine 


Fig. 1 Fig. 2 


Zum Verhalten des Torsions- Die Knicklasten des durch ein semitangentiales Moment 
momentes beim Ausknicken. und eine axiale Druckkraft beanspruchten Stabes fiir ver- 


schiedene Arten der Lagerung («% = Biegesteifigkeit, 
1 = Lange des Stabes). 


Belastung einen Punkt innerhalb der zu seiner Lagerung gehérenden kritischen 
Kurve ergibt. Die Gleichungen der Kurven sind fiir die Lagerungen I, IV, V | 
transzendent, lassen sich aber durch einfache Polynome gut annahern. Da in den 
Fallen I und IT der belastende Momentvektor seine Richtung nicht andern kann, 
sind die zugehérigen Kurven allgemein giiltig. 

Ein Bericht tiber diese Untersuchungen und eine Zusammenstellung der Re- 
sultate sollen bald veréffentlicht werden. 


Errata 


Die Knicklast des einseitig eingespannten, tangential gedriickten 
Stabes. Von Max Breck (ZAMP 3, Fasc. 3, 225 [1952]). 


Auf Seite 228, zweite Zeile, ist der Satz «Im Falle einer Doppelwurzel kann 
y die Form 
y(¥, t) = Y(x) tsin(o 2) 
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annehmen » zu ersetzen durch: «Tritt fiir P = P, eine Doppelwurzel auf, so treten 
bei einer beliebig kleinen Erhéhung von P iiber #, hinaus an Stelle der zwei 
verschwundenen reellen zwei komplexe Wurzeln fiir w?2 1), und y kann die Form 


y(#, t) = Y(x) e®”* sin (o’ 2) (w” > 0) 


annehmen.» — Auf die iibrigen Teile der Arbeit und insbesondere auf die Schluss- 
folgerungen wirkt sich diese Anderung nicht aus?). 


Uber die Koppelkurven des Zwillingkurbeltriebes. Von W. MEYER zuR 
CAPELLEN (ZAMP 2, Fasc. 3, 189-207 [1951)). 


Seite 189: Dass gewisse Koppelkurven des Antiparallelkurbeltriebes als Um- 
riss von Tragfliigelprofilen benutzt werden kénnen, wurde nicht von PIERCY usw. 
(Anmerkung 1, Seite 189) angegeben, sondern wohl zuerst von CLEFF und WEL- 
BOURNE (beide England) ausgesprochen. Den Beweis fiir die Identitat dieser 

_ Koppelkurven mit den zu den Kegelschnitten inversen Kurven gab O. S. HEckx, 
Z. angew. Math. Mech. 30, 239 (1950). 
Seite 194, Zeile 18 von unten: Statt # sint muss es (c sint)/2 heissen. 
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Naturforschung und Medizin in Deutschland 1939 bis 1946 (fiir den 
Buchhandel bestiminte Ausgabe der Fiat Review of Geyman Science) (Dieterich- 
sche Verlagsbuchhandlung, Wiesbaden). Preis je Band DM. 10.—. 


Band 8 und 9: Physik der festen Korper. Herausgegeben von G. Joos (Wies- 
baden 1947 bzw. 1948), 228 und 235 S. 

In diesen zwei Banden berichten achtundzwanzig durchwegs wohlausgewiesene 
Autoren in zusammenfassenden Referaten iiber Fortschritt und Stand der For- 
schung in den verschiedenen Zweigen der Festkérperphysik. Der Inhalt umfasst 
in grossen Ziigen: Struktur, mechanische Eigenschaften, Thermodynamik, 
magnetische Eigenschaften, elektrische Eigenschaften und Optik der festen 
Korper. Jedes dieser Kapitel ist durch mehrere griindliche und konzentriert ge- 
haltene Aufsatze belegt. ; 

Wenn auch in einzelnen Abschnitten die Darstellung den Zusammenhang mit 
den Forschungen in andern Landern zu sehr vernachlassigt, SO erfullt doch 
das Buch seinen Hauptzweck aufs beste: Es bringt eine griindliche Ubersicht 
iiber das wahrend der Kriegszeit nur liickenhaft zu uns gelangte, umfangreiche 
deutsche Schrifttum iiber Festkérperphysik. H, Labhart 


Bd. 13: Kernphysik und kosmische Strahlen, Teil. I. Von W. BorHEe und 
S. FLUGGE (Wiesbaden 1948), 230 Seiten. 

In den beiden ersten Abschnitten des Buches wird der Problemkreis der 
Elementarteilchen und der Felder, der vor allem fiir die Elektrodynamik und 


1) Vergleiche Ernst LINDELOF, Le calcul des résidus (Gauthier: Villars, Paris 1905), 3 ae 
oder Epovarp Goursat, Cours @analyse mathématique, Bd. 2, Kapitel XV Il: Théoréme de Weier- 
strass. - Der Ausdruck « Doppelwurzel» wird im Hinblick aui Figur 4 een acct 

2) Die Aufdeckung dieses Fehlers verdanke ich Herrn Prof. Dr. R. von Mises. 
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die kosmische Strahlung von fundamentaler Bedeutung ist, sowohl vom theore-_ 
tischen als auch vom experimentellen Standpunkt aus behandelt. Die restlichen 
Kapitel befassen sich mit der Physik der Atomkerne, der Kernumwandlungs-— 
prozesse und der Spaltung der schwersten Kerne. Die Abhandlungen geben den — 
deutschen Beitrag zu den aktuellen Problemen der Physik wieder; dabei mussten 
zufolge des fehlenden Kontaktes zahlreiche auslandische, vor allem amerikanische 
Arbeiten unberiicksichtigt bleiben. Trotzdem sind die Berichte stellenweise fiir 
den Fachmann fesselnd geschrieben. H. Waffler 


Einfiihrung in die Differentialgeometrie. Von W. BLASCHKE (Springer-_ 
Verlag, Berlin 1950). 146 S., 57 Abb.; DM 16.-. { 


Das neue Lehrbuch von WILHELM BLASCHKE iiber Differentialgeometrie ver-_ 
mittelt eine ausgezeichnete, moderne Einfiihrung in dieses interessante Gebiet 
der Geometrie, dessen Anfange auf C. F. Gauss und G. MonGcE zurtickgehen. 
Eingangs sind auf knappem Raume die Hilfsmittel aus der linearen Algebra — 
zusammengestellt, anschliessend folgt die Behandlung der Streifen und Linien 
mit Hilfe des begleitenden Dreibeins (repére mobile). Der dritte Teil bringt eine © 
Einfiihrung in den Cartanschen Kalkiil der alternierenden Differentialformen 
(Pfaffsche Formen), der fiir die folgenden Entwicklungen grundlegend ist. Die © 
innere Flachenlehre, welche die bei einer Verbiegung invarianten Eigenschaften 
der Flache umfasst, ist im nachsten Kapitel dargestellt, an das die Untersuchung ~ 
der geodatischen Linien ankniipft. Alsdann folgt die 4ussere Flachenlehre (bewe- — 
gungsinvariante Eigenschaften); und schliesslich die Behandlung der Minimal- 
flachen. ' 

Jedem Kapitel sind eine Reihe von zum Teil neuen Aufgaben, die einen © 
wesentlichen und wertvollen Bestandteil des Buches bilden, sowie Erginzungen — 
zur Theorie beigefiigt. Das bemerkenswerte Werk wird wirksam belebt und abge- 
rundet durch zahlreiche in den Text eingestreute historische Bemerkungen und 
Hinweise auf aktuelle Probleme der Forschung. Die Darstellung, die sich auf das 
Wesentliche beschrankt, ist durch die konsequente Verwendung des Cartanschen — 
Kalkiils sehr elegant und kurz, 

Diese meisterhafte Einfiihrung in die Differentialgeometrie kann allen Inter- — 
essenten, zu denen neben den Mathematikern auch die Geodaten, Astronomen 
und Physiker zu zahlen sind, sehr empfohlen werden. E. Roth-Desmeules 


—e « 


Linear Computations. By P. S. DwyER (John Wiley & Sons, New York 
1951), 344 pp.; $6.-. 


Das Buch will Statistiker mit den modernen Methoden zum Auflésen linearer 
Gleichungssysteme bekanntmachen. Ein Einfiihrungskapitel hat das Rechnen 
mit Annaherungswerten (Fehlerrechnung) und andere grundsatzliche Fragen des 
numerischen Rechnens zum Gegenstand. Darauf folgen langere Ausfiihrungen 
iiber die Auflésung von linearen Gleichungssystemen und die Inversion von 
Matrizen, beginnend mit den zum Teil veralteten Verfahren wie «Multiplikation-— 
Subtraktion» bis zu den modernen Matrizenverfahren. Dabei sind zahllose nume- 
rische Beispiele in den Text eingeflochten, die besonders fiir den Anfanger sehr 
instruktiv sind. Ferner enthalt das Buch einige Abschnitte tiber die Berechnung 
von Determinanten, Eigenwerten sowie tiber Matrizenrechnung im allgemeinen. 
Leider wird aber auf iterative Methoden wie die Relaxation tiberhaupt nicht ein- 
gegangen. Den Abschluss des Buches bildet ein Kapitel mit zahlreichen Anwen- 
dungsbeispielen aus der Statistik. HA. Rutishauser 
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